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A. Introduction 


This document is intended to provide the specifications and requirements for a 
flight control system design challenge. The response to the challenge will involve 
documenting whether the particular design has met the stated requirements through 
analysis and computer simulation. The response should be written in the general format 
of a technical publication with corresponding length limits, e.g., an approximate 
maximum length of 45 units, with each full-size figure and double-spaced typewritten 
page constituting one unit. 


B. Challenge Goal and Vehicle Description 



The underlying goal of the challenge is the design and computer simulation of a 
high-performance piloted aircraft in a series of different flight conditions. Of particular 
importance to the design will be the ability of the aircraft to sustain significant system 
failures/damage, here to be represented by variations in control surface actuator 
characteristics and vehicle dynamics. The vehicle that is the subject of the challenge is 
the Innovative Control Effector (ICE) vehicle shown in Fig. 1.1 . A detailed description of 
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Figure 1. 1 The ICE vehicle 
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and the derivation of a linear model is given in the Appendix. 

C. Design Challenge 

1 . Desired Control System Response Types 

The control system response types are straightforward: For longitudinal control, a 
pitch-rate command system should be created. There is no requirement for attitude hold 
in this system. For lateral control, roll-rate and sideslip-command systems should be 
created, again with no requirement for attitude hold in the roll rate system. Commands to 
the pitch-rate and roll-rate systems are assumed to be generated by cockpit inceptor 
commands, e.g., longitudinal and lateral control column inputs. No pilot-generated 
sideslip commands will be in evidence, i.e., the system should maintain zero sideslip with 
the pilot’s “feet on the floor”. The vehicle models that will be provided have been 
linearized about a series of different flight conditions. The maneuvers to be described 
will be transient in nature and no thrust changes from trim values will be required. 

The desired characteristics of the rate command systems will be predicated upon 
meeting handling qualities requirements in the pitch and roll axes to be specified in 
Section I.C.5.2. However desired frequency-domain characteristics can be specified as: 

_9 100 

0 C s(s + 5)(s + 20) 

w — all for 0.5 < to < 50 rad / sec ( 1 ) 

<t> Sc s(s + 5)(s + 20) 

(3 100 

p c ~ (s + 5)(s + 20) 

where 0 and 0 C represent the Euler pitch attitude and attitude commands, and <t> s and 4> Sc 
represent roll attitude and attitude commands about the x-stability axis, i.e. describing 
rolling motion about the vehicle velocity vector in equilibrium flight. 

Finally, in keeping with the design philosophy of the ICE vehicle, the controller 
must utilize all the available control surface effectors. 

2. Controller Simulink Implementation 

The controller must be implemented as a discrete device in the Simulink computer 
simulation that will form an important part of the Design Challenge response. A frame 
rate of 80 Hz may be assumed in the discretization. 
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3. 


Pilot Models 


A pair of analytical pilot models will be specified as part of the computer simulation 
of the pilot/vehicle system. Figure 1.2 is the Simulink® diagrams for the pilot models. 
No pilot model for the beta loop is necessary because of the assumption of “feet-on-the- 
floor” activity by the pilot. The pilot models of Fig. 2 have been created assuming the 
dynamics of Eq. 1 are in evidence in the 0 and <f> s loops. The gains of 2.67 on visual error 
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Figure 1.2 Pilot models 


in each model will produce crossover frequencies of 1 .0 rad/sec in each of the loops, a 
value representative of moderate pilot control activity. Figure 1.3 is the Bode plot of the 
open-loop transfer function of the pilot/vehicle system for both the 0 and (f>-loops when 
the dynamics of Eq. 1 are in evidence. 
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Figure 1.3 Bode diagram of nominal pilot/vehicle open-loop system 
3. Sensor Models 

It will be assumed that any and all of the output variables specified in the 
Appendix for the linear models to be discussed will be available for feedback in the 
control system design. These include v, a w qt,, 0, Pw, ps, r s , <|>, a Xcg , a ycgj an d a ncg No 
sensor dynamics will be modeled with the exception of a 0.02 sec time delay assumed for 
each sensed variable. Additive sensor noise must be included in each measured variable 
as shown below in Fig. 1.4. 
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Figure 1.4 Sensor noise models 


The noise scale factor is given as K = 0.15625 for all sensed variable except for 0 
and <b, for which K = 0, and for a x , a v , and a n for which K = 0.015625. The 

T 5 3 A Cg 7 /Qg Cg 

variables for 0 and <j> should not be differentiated in the control law to be determined. 
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With K = 0.15625, the root mean square (RMS) value of the output of the filter in Fig. 1.4 
will be 0.25. Thus, for example, the RMS noise on a w would be 0.25 deg, etc. 

4. Linear Models - Flight Conditions 

Linear models for the ICE vehicle are enclosed in the Appendix for the following 
four flight conditions: 

Level Flight: Mach No. = 0.3 Altitude = 15,000 ft 

Mach No. - 0.6 Altitude = 25,000 ft 
Mach No. = 0.9 Altitude = 35,000 ft 

Steady Turning Mach No. = 0.6 Altitude = 25,000 ft 
Flight 

The latter condition involves a trim bank angle of 60 deg. 

5. Piloting Tasks 

The piloting tasks will consist of the pilot following a series of simultaneous pitch 
and roll attitude commands for a period of 50 sec in each flight condition. The 
commands consist of a series of filtered 5 sec pulses, alternating in sign as shown in Fig. 
1.5. The amplitude of the pulses depends upon the flight condition and the command, i.e., 
whether it is a pitch or roll command. As will be described in Section I.C.8, in 
comparing the command input to the corresponding pilot/vehicle response, a 1.0 sec 
delay can be added to the recorded command signal. This is only for the purposes of 
performance assessment. The delay is not included in the input to the pilot/vehicle 
system in the Simulink® simulation. Figure 1.6 shows the command generator. 
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Figure 1.6 Command signal generation 
The amplitudes of the filtered pulses are given as 

Level Flight: Mach No. = 0.3 Altitude = 15,000 ft Gamp = 5 deg; <f>am P = 30 deg 

Mach No. = 0.6 Altitude = 25,000 ft Gamp = 5 deg; <t»am P = 30 deg 

Mach No. = 0.9 Altitude = 35,000 ft Gam P = 0deg; 4» amp = 5 deg 

Steady Mach No. = 0.6 Altitude = 25,000 ft Gam P = 0 deg; 4>amp = 5 deg 

Turning 

Flight 

6. Performance/Handling Qualities Specifications 

a. Desired Handling Qualities 

Handling qualities specifications are to be analytically determined using 
bandwidth/phase delay as applied to the pitch-rate and roll-rate systems. Figure 1.7 is to 
be used to define bandwidth and phase delay, and Fig. 1.8 is to be used to assign handling 
qualities levels. NB: The transfer functions to be used in applying the measures of Fig. 
1.7 are G/0 Cps and <t>/<|>c FS where G Cps and <|> CFS include the dynamics of the force/feel system 

given by 

625 

force /feel dynamics = — (2) 

s 2 + 35s + 625 

Note that the units on the force/feel system dynamics are of no importance in this 
analysis. 
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Figure 1.7 Determining bandwidth for handling qualities determination 
Bode diagram is for either 0/0c FS or <|>/<t>c FS 


Phase delay t p is determined from Fig. 1.7 as: 

AO 

Tp _ (57.3)(2 x co ]g0 ) 


( 3 ) 



Figure 1.8 Handling qualities level - bandwidth/phase delay mapping for pitch axis; for 
roll axis, the roll bandwidth determined from Fig. 7 shall be at least 1 rad/sec for Level 1. 
Phase delay shall be no more that 0.14 sec for Level 1 and 0.2 sec for Level 2 
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b. Desired Tracking Performance 

The pilot/vehicle performance specifications in each of the flight conditions are: 

pitch-attitude tracking: 

Gamp = 5 deg or 0 deg 

Desired: no sustained oscillations; maximum error of ± 2 deg 

Adequate: no sustained oscillations: maximum error of ± 3 deg 

roll-attitude tracking: 

<j>amo = 30 deg 

Desired: no sustained oscillations; maximum error of ±10 deg 

Adequate: no sustained oscillations; maximum error of ±15 deg 

^amo ~ 5 deg 

Desired: no sustained oscillations; maximum error of ± 2 deg 

Adequate: no sustained oscillations; maximum error of ±3 deg 

sideslip tracking: 

All attitude command amplitudes 

Desired: maximum excursions of ± 2 deg 

Adequate: maximum excursions of ± 5 deg 


c. Desired Stability Margins 

Stability margins in the control system design are to be determined in classical 
single-loop fashion by sequentially cutting each control loop before each actuator as 
indicated in Fig. 1.9. 



Figure 1.9 Determining loop break point 
X 

From the Bode diagram of -^-(jto) the following stability margins should be in 
evidence: Gain Margin > 6 dB; Phase Margin > 30 deg 
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NB: It is obvious that this stability analysis is not as rigorous as one obtained by 

inserting a perturbation matrix of the form P = diag^^ - * 1 ,...,K n e~*" ) before the 

actuators and assessing closed-loop stability when Kj and <(>j are varied within some 
desired region in the gain and phase parameter space. The single-loop approach was 
adopted for the sake of simplicity here given the number of actuators involved in the 
vehicle. 

d. Structural Coupling Considerations 

Although no elastic modes have been included in the linearized ICE models, to 
minimize excitation of structural modes the following requirements on closed-loop 
transfer functions should be met for <o < 25 rad/sec: 
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( 4 ) 


e. Control Activity 

In each of the flight conditions for the tracking tasks to be considered, control 
activity should meet the following criterion: 

The half-power frequency in the command signal to each actuator should be less 
than or equal to the actuator bandwidth (here assumed equivalent to the undamped 
natural frequency of the actuator. 

In equation form, this requirement can be stated as 

(C0 05 ), <(()>„), (5) 

where (coo s)i is obtained as 

(<°o s)i 

0.5= J0 5 8 (to)d(o (6) 

0 “ c ' 

with 8 Cj denoting the command to the i 111 actuator and <D xx (g>) denoting the power spectral 

density of the signal x(t). Finally, (co n )j denotes the undamped natural frequency of the ith 
actuator. The requirement of Eq. 5 places a limit on the control activity required by the 
control system and the sensor noise transmission that occurs as a result of the design. 
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f. Scheduling Considerations 

In the analysis and computer simulation of the pilot/vehicle system, the flight 
control design may be tuned to each flight condition being considered. The one 
exception to this is the turning flight condition, where the level flight dynamics must be 
considered in the control system design. No scheduling algorithms need be developed. 
However, a discussion of the approach to scheduling should be included. Finally, to 
demonstrate the robustness of the design approach, the pilot/vehicle system tasks should 
be performed at each of the four flight conditions using the control system designed for 
the Mach 0.6 Altitude 25,000 ft flight condition. No handling qualities evaluations need 
be included, however, tracking task performance should be ascertained using the 
requirements of Section l.C.7.b. 

7. Failure Modes 


After computer simulation of the nominal system in the specified flight 
conditions, the following system failure(s) should be evaluated via Simulink® simulation. 
These failures should be evaluated at each flight condition with the tracking tasks 
specified in Section I.C.5. 


Failure No. 


Failure Description 


1 left elevon actuator amplitude limits reduced to ±1 5 deg 
rate limits reduced to ±1 5 deg/sec 

2 symmetric pitch flap actuator with hard-over of +5 deg 

3 left spoiler actuator with rate limits reduced to ±10 deg/sec 

4 left leading edge outboard flap with hard-over of +5 deg 

5 pitch nozzle actuator with rate limits reduced to ±10 
deg/sec 

6 elements of linear plant A and B matrices each changed by 
±20 %. This change should not include those elements of 
the A matrix that describe kinematic relationships. To 
ensure uniformity, the sign of the 20% variations will be 
dependent upon the column of the A and B matrices, with 
even- numbered columns receiving a ±20% change and 
odd- numbered columns receiving a - 20% change. 


It is desired that the pilot/vehicle system be able to tolerate all six of these failures 
simultaneously. In the Simulink® simulation, the failures should be introduced at t = 20 
sec. Section I.C.8 briefly defines an M-file that will automate the failures. 
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In the event that the simultaneous failures prove too severe, a series of Simulink 
simulations can be undertaken in which the failures are added incrementally. If this 
procedure is followed, the order in which the failures are accumulated should follow the 
failure numbers given above. That is, first, No. 1, then Nos. 1 and 2, etc. 

8. MATLAB® Models and M-file for Simulink® Failure Simulation 


A MATLAB® M-file called “IceFunctionsDC.m” is provided that automates the 
failures just described. The file can be used with a Simulink® simulation module called 
ICEdc.mdl containing the state-space vehicle model, actuator models, input commands, 
pilot models and reference models. Figure 1. 10 shows ICEdc.mdl. 
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Figure 1. 10 The ICEdc.mdl in Simulink' 




By left clicking on any of the blocks above, the underlying model structures can be 
viewed. By clicking on the shadowed box containing “DC = Wings Level...” the menu 
shown in Fig. 1.1 1 appears. 



Figure. 1.1 1 ICEdc.mdl menu 
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The menu of Fig. 1.1 1 allows the user to select flight conditions, and “design” conditions. 
If desired, the user can load their control designs into MATLAB, using the “Other Design 
Conditions” label. Likewise, the Flight Conditions can be loaded into MATLAB using 
the “Model File Name.” For example, in Fig. 1.11, the “Design Condition” is Wings 
Level, Mach No. = 0.3 and Altitude = 15,000 ft, i.e. this is the flight condition for which 
the control system has been designed or “tuned”. The flight condition selected is also 
Wings Level, Mach No = 0.3 and Altitude = 15,000 ft. The octagonal elements on the 
right hand side of the model will terminate the Simulink® simulation when the associated 
output variables exceed the values within the elements, e.g., here 50, indicating 
instability. The Model Reference is also included. Depending upon the particular design 
philosophy adopted by the user, this block may or may not be necessary, i.e. the reference 
model dynamics may be implicit in the design. 

Returning to Fig. 1.10, one initializes the system by left clicking on the shaded box 
labeled INITialize Nominal Model. This sets the actuators to nominal conditions and the 
vehicle model to the flight condition specified in the menu of Fig. 1.1 1. By left clicking 
on the Fail Trigger box, one can select a failure time, now set at a default value of 20 sec. 
Alternatively, the user can manually fail the system by left clicking on the FAIL model 
box any time during the simulation run. The failure conditions (including the ±20% 
variation in appropriate elements of the A and B matrices) must be specified by the user 
through the m-file IceFunctionsDC.m. listed in Appendix I.. 

9.0 Design Challenge Summary 

The Design Challenge can be summarized as follows: 

Vehicle: ICE aircraft, linearized about four flight conditions. 

Controller: Designed by participant. May be tuned to each of the four flight 

conditions, with the exception that for the turning flight condition, the 
controller is tuned to the wings-level condition. Each of the control 
effectors must be utilized in the design. Controller is to be discretized for 
Simulink® computer simulation. 

Handling Qualities: Ascertained through bandwidth/phase delay for each of the four 

flight conditions, for nominal vehicle only (no failures). 

Tracking Performance: Ascertained through Simulink® computer simulation for each 

flight condition. Should include nominal and failed cases. 
Failure should occur 20 seconds into a 50 second run. 

Stability Margins: Ascertained through linear analysis for each of the four flight 

conditions, for nominal vehicle only (no failures). 
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Structural Coupling: Ascertained through linear analysis for each of the four flight 

conditions, for nominal vehicle only (no failures). 

Control Activity: Ascertained through Simulink® computer simulation for each 

flight condition, for nominal vehicle only (no failures). 

Scheduling Considerations: Scheduling need not be implemented, but means for 

accommodating scheduling (if necessary) should be 
discussed. 

Off-nominal Flight Cond’s: Ascertained through Simulink computer simulation for 

each of the four flight conditions. Controller designed for 
Mach No. = 0.6, Altitude = 25,000 ft should be employed 
in remaining three flight conditions. Tracking 
performance should be assessed. 

D. Appendix I Listing for IceFunctionsDC.m 


function [errFlag,tstep,sYsDC / SYsFC / sYsFail,icF[:,icDC] = IceFunctionsDC (inVar ) ; 
k 

VThis m-file contains all required functions for the ICE Simulink model 
* 

^Inputs: 

k inVar - input variable, may be a string or numeric vector. The following cases are defined: 
** 'Init 1 - Update model (3tate space and actuators) with nominal system 

k ’Fail' - Update model (state space and actuators) with failed system 

k [0 failFlagJ - fail trigger input vector 

% [ia xl x2 ...] - call various U3er defined functions with arguments xl x2 ... 


k Outputs: 

% errFlag - error flag, G=no error, l=error 
k tstep - simulation time step 

^ SysFC - system state space matrices for current flight condition 
k SysDC - Design condition state space matrices 
k SysFail - failed system state space matrices 

k uDC,xDC,yDC - Initial conditions and inputs for Design condition 
k SIMULINK BLOCKS - various blocks updated for each different case 

currentSys = gcs; ^parent name of current SIMULINK object 

slashSpot - findstr (currentSys, 1 / 1 ) ; *find fwd slashes in current sys object- 
if size (slashSpot, 2) > 0 

fileName = currentSys (1: slashSpot (1) -1) ; ^delete everything to left of last slash 
else 

fileName - currentSys; 

end 
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k Determine which function to switch to 
if ischaE (inVar) k±3 the input a character or a number 
if strcmp (inVar , 1 Ini t’ ) 

myCase = 1; * Initialization of Nominal model 

elseif strcmp (inVar Fail 1 ) 

myCase = 2; ^Initialization of Failed model 

else 


myCase » 0; 

end 

elseif inVar (1)==0 
if inVar (2) ==1 
myCase = 2; 
else 

myCase =3; 
end 

elseif inVar (1)==1 
myCase = 4; 
elseif inVar (l)-=2 
myCase = 5; 
elseif inVar ( 1 ) = =3 
myCase = 6; 
else 

myCase = 0; 


% Invalid input string 

ki.3 input a number ?H Check fail trigger 
*IF Fai lFlag==true 
^Fail the vehicle 

*sDo nothing, no error flag returned 

HVo Something per your design. 

HDo Somethin per your design. 

^Po something per your design. 
^Invalid input 


end 

errFlag = 1; 


k Output. Variable D ef inti 0213 

tstep = 0.0005; % simulation time step 

k Initialise system state space matrices for designed condition 
k This is only included for possible aid in simulating flight evelope testing. 

% It does nothing in the given SIHULINK file! 

DC = get_j)aram (str cat ( filename, 1 /Flight Condition’ ), 'DC ') ; *Read desired flight condition 
dsgnFile=get_param ( street (fileName, '/Flight Condition' ),' dsgnFile ' ) ; %Read posible input fil 
if strcmp (DC, ' Other ...') *=1; *if Dropdown says "other” right dsgnFile to DC 
DC=dsgnFile ; 
end %if 

k Initialize system state space matrices for flight condition 

k This writes to the block FSAVmodel whatever model is in the Flight dropdown or if "other 
k is selected then it reads the input file 
FC « get_par am (str cat ( f ileName, 1 /Flight Condition' ),' FC 1 ) ; 
modFile=get_param(3trcat(fileName, ' /Flight Condition' ) , 'ModelFile ‘ ) ; 
if strcmp (FC, ' Other ... 1 ) ==1; 

FC=*modFile ; 
end kit 

k Call function to write FC to model, DC only writes DC to sysDC nothing elese 
[sysFC,sysDC , ays Fail ,icFC, icDC]=getICEsystem(DC,FC) ; 

switch myCase 
case 0 

errFlag = 1; 
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NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL NOMINAL 


case 1 

^ Initialize Nominal system and actuators in Simulink 
% Inputs: 

^ A,B,C,D - nominal state 3pace matrices 
k fileName - current SIMULINK file name 
% Outputs: 

* SIMULINK OBJECTS - update state 3pace system block and all actuator components 
*Left Elevon (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 


delay 

= 

0.00001; 

^Transport delay (sec) 

k 

= 

1; 

<gain 

w 

= 

63.2455532; 

%freqency (rad/s) 

damp 

» 

1.106797181; 

i- damping ratio 

dMaxHi 

= 

30; 

imax position, upper limit 

dMaxLo 

= 

-30; 

-;max position, lower limit 

ddotHax 

= 

150; 

imax rate (+/- unit/sec), enter positive number 

hardOver 

= 

0; 

%hard-over switch (0=no hard-over, l=hard-over) 

jamPos 

= 

0; 

^hard-over actuator position (unit) 

setActuators (fileName, ‘03 

‘ , [delay;k ;w; damp; dMaxHi ; dMaxLo; ddotHax ;hard0ver; jamPos] ) ; 


%Right Elevon (2nd Order, with limits); G(s) » k * w A 2/(3 A 2 + 2*damp*w*s + w A 2) 
delay * 0.00001; ^Transport delay (sec) 

k =1; isgain 

w - 63.2455532; ^freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi - 30; *max position, upper limit 

dMaxLo * -30; ^max position, lower limit 

ddotHax = 150; ^max rate (4-/- unit/sec) , enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos = 0; ^hard-over actuator position (unit) 

setActuators (fileName, 1 13 ' , [delay;k; w; damp; dMaxHi; dMaxLo; ddotHax ;hardOver; jamPos] ) ; 

^Symmetric Pitch Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay * 0.00001; ^Transport delay (sec) 

k =1; ^gain 

w = 63.2455532; %£reqency (rad/s) 

damp - 1.106797181; ^damping ratio 

dMaxHi = 30; ^max position, upper limit 

dMaxLo = -30; %max position, lower limit 

ddotHax * 50; ^max rate (+/- unit/sec), enter positive number 

hardOver = 0; %hard-over switch (0=no hard-over, l=hard-over) 

jamPos « 0; ^hard-over actuator position (unit) 

setActuators (fileName, ’4* , [delay;k;w;damp;dMaxHi;dMaxLo;ddotMax;hardOver; jamPos] ) ; 
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^Left All Moving Tip (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2 1r damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; *gain 

w = 63.2455532; %fregency (rad/s) 

damp = 1.106797181; ^damping ratio 

dHaxHi = 60; *max position, upper limit 

dMaxLo = 0; ^max position, lower limit 

ddotMax = 150; %max rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos ■ 0; ^.hard-over actuator position (unit) 

setActuators (f ileName, ' 05 ' , [delay;k;w;damp;dMaxHi;dMaxLo;ddotMax;hardOver; jamPos] ) ; 

%Right All Moving Tip (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*3 4- w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k * 1; %gain 

w = 63.2455532; ^freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi = 60; %max position, upper limit 

dMaxLo = 0; ^max position, lower limit 

ddotMax = 150; %max rate (+/- unit/sec), enter positive number 

hardOver = 0; Tfhard-over switch (0=no hard-over, l=hard-over) 

jamPos = 0; %hard-over actuator position (unit) 

setActuators (f ileName, 1 15 1 , [delay;k;w;damp;dMaxHi;dMaxLo;ddotMax;hardOver; jamPos] ) ; 


^Left Spoiler (2nd Order, with limits); G ( 3 ) * k * w A 2/(s A 2 + 2*damp*w*3 + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; ^gain 

w = 63.2455532; *freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi = 60; ^max position, upper limit 

dMaxLo « 0; %max position, lower limit 

ddotMax » 150; %max rate (+/- unit/sec) , enter positive number 

hardOver = 0; Vhard-over switch (D=no hard-over, l=hard-over) 

jamPos = 0; \hard-over actuator position (unit) 

setActuators (f ileName, 1 09 1 , [ delay; k ;w; damp; dHaxHi ; dMaxLo; ddotMax ;hard0ver; jamPos] ) ; 

s ;Right Spoiler (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k = 1 ; ^gain 

w = 63.2455532; *freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi = 60; hmox position, upper limit 

dMaxLo =* 0; %max position, lower limit 

ddotMax = 150; ^max rate (4/- unit/sec) , enter positive number 

hardOver * 0; %hard-over switch (0=no hard-over, l=hard-over) 

jamPos = 0; *hard-over actuator position (unit) 

setActuators (f ileName, 1 19 1 , [del ay;k;w; damp; dHaxHi ;dHaxLo; ddotMax ;hard0ver; jamPos ] ) ; 
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%Left Leading Edge Outboard Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp 1c w*s + w A 2) 
delay - 0.00001; ^Transport delay (sec) 

k * 1; %gain 

¥ = 42.22321636; ^freqency (rad/s) 

damp ® 1.3952987; ^damping ratio 

dMaxHi * 40; %max position, upper limit 

dMaxLo = ^40; ^max position, lower limit 

ddotMax = 40; *max rate (+/- unit/3ec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos * 0; ^hard-over actuator position (unit) 

setActuators (fileName, 1 02 ' , [ del ay;k;w; damp; dMaxHi; dMaxLo; ddotMax ;hard0ver; jamPos] ) ; 

%Right Leading Edge Outboard Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay * 0.00001; ^Transport delay (sec) 

k =1; *gain 

w * 42.22324636; *freqency (rad/s) 

damp = 1.3952987; ^damping ratio 

dMaxHi = 40; ^max position, upper limit 

dMaxLo = -40; *max position, lower limit 

ddotMax = 40; ^max rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l^hard-over) 

jamPos = 0; Tihard-over actuator position (unit) 

setActuators (fileName, 1 12' , (delay;k;w;damp;dMaxHi;dMaxLo;ddotMax;hardOver; jamPos] ) ; 


^Pitch Nozzle (2nd Order, with limits); G(s) * k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay ■ 0.00001; ^Transport delay (sec) 

k - 1; ^gain 

w * 39.18905204; ^freqency (rad/s) 

damp ■ 1.001376608; ^damping ratio 

dMaxHi = 15; ^max position, upper limit 

dMaxLo - -15; ^max position, lower limit 

ddotMax ■ 60; %max rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) NOTE: INTEGER 

jamPos = 0; ^hard-over actuator position (unit) 

setActuators (fileName, ' 10 ' , [ delay;k;w; damp; dMaxHi; dMaxLo ;ddotMax;hard0ver; jamPos] ) ; 

Warn Nozzle (2nd Order, with limits); G(s) - k * w A 2/(3 A 2 + 2*damp*w*s + u A 2) 
delay * 0.00001; ^Transport delay (sec) 

k = 1; ^gain 

w * 39.18905204; %f regency (rad/s) 

damp = 1.001376608; ^damping ratio 

dMaxHi = 15; %max position, upper limit 

dMaxLo ® -15; ^max position, lower limit 

ddotMax =* 60; ^max rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) NOTE: INTEGER 

jamPos = 0; ^hard-over actuator position (unit) 

setActuators (fileName, ' 20 ' , [ delay ;k;w; damp; dMaxHi; dMaxLo; ddotMax; hardOver; jamPos] ) ; 
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set_param(strcat (f ileName, 1 ./Actuators 1 ) , 'EackgroundColor ' , ‘white 1 ) 

^Initialize the 3tate space model and Input amplitudes 

setjparam (str cat (f ileName, ' /ICEmdl ' ) , 1 A ' ,»at2str (sysFC. a, 6) , 'B 1 ,mat2str (sysFC.b,6) , 1 C ' , 
mat2str ( sysFC. c, 6) , 'D ’ ,mat2str (sysFC. d,6) , 'XO ' ,mat2str (icFC.xO,6) ) ; 
set_param ( street (f ileName, 1 /ICEmdl ' ) , 'EackgroundColor 1 , 'white ' ) ; 
set_param(strcat (f ileName, ' /InputCmds/thetaamp 1 ) , 'Value' ,num2str (icFC. thetaamp) ) ; 
3et_param (streat (f ileName, ' / Input Cmds /phi amp 1 ) , 'Value 1 ,num2str (icFC.phiamp) ) ; 

disp ( ' IceFunctionsDC» Initialization of model completed...') 

disp(' Output 3 [ errFlag, tstep,sysDC, sysFC, sysF&il,icFC,icDC] written. ') 

errFlag = 0; 


* FAILED FAILED FAILED FAILED FAILED FAILED FAILED FAILED FAILED FAILED FAILED FAILED 


case 2 

k Initialize FAILED 3y3tem and actuators in Simulink 
k Inputs: 

k Af ail ,Bf ail, C fail ,D fail - FAILED state space matrices 

k fileNaiae - current SIMULINK file name 
k Outputs: 

k SIMULINK OBJECTS - update state space system block and all actuator components 


^FAILED Left Elevon (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*u*s + w A 2) 


delay 


0.00001; 

^Transport delay (sec) 

k 

= 

1; 

*gain 

w 

= 

63.2455532; 

kt regency { rad/s) 

damp 

= 

1.106797181; 

% damp mg ratio 

dMaxHi 

= 

30; 

%max position (unit), upper limit (Nom:30) 

dMaxLo 

= 

-30; 

*max position (unit) , lower limit (Nom:3D) 

ddotMax 

= 

15; 

*max rate {+/- unit/sec), enter positive number (Nom: 150) 

hardOver 

= 

0; 

^hard-over switch (0=no hard-over, l=hard-over) 

jamPos 

= 

0; 

^hard-over actuator position (unit) 

setActuators (f ileName, '03' 

' , [delay;k;w;damp;dMaxHi;dMaxLo;ddotMax;hardOver; jamPos] ) ; 


^FAILED Right Elevon (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; ^gain 

w = 63.2455532; %freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi = 30; ^max position, upper limit 

dMaxLo * -30; ^max position, lower limit 

ddotMax = 150; %max rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos * 0; ^hard-over actuator position (unit) 

setActuators (f ileName, ’ 13 ' , [delay;k;w;damp;dHaxHi;dF!axLo;ddotJ!ax;hardOver ; jamPos]) ; 
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^FAILED Symmetric Pitch Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*3 + u A 2) 

delay » O.OOOOl; ^Transport delay (sec) 

k = 1 ; %gain 

v = 63.2455532/6.32455532; %freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dHaxHi = 30; ^max position, upper limit 

dHaxLo = -30; ^max position, lower limit 

ddotHax » 50; ^max rate (+/- unit/sec), enter positive number 

hardOver = 1; ^hard-over switch (0-no hard-over, l=hard-over) 

jamPos = 5; ^hard-over actuator position (unit) 

3 etActuator s ( f ileName, 1 4' , [delay;k;w;damp;dMaxHi;dMaxLo;ddotHax;hardOver; jamPos] ) ; 

* FAILED Lett All Moving Tip (2nd Order, with limits); 0(s) = k * w A 2/(s A 2 + 2*damp*w*s + u A 2) 
delay = 0.00001; ^Transport delay (sec) 

k * 1; *gain 

w = 63.2455532; ^freqency (rad/s) 

damp ■ 1.106797181; ^damping ratio 

dHaxHi = 60; *max position, upper limit 

dHaxLo = 0; ^wax position, lower limit 

ddotHax = 150; smax rate (+/- unit/sec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, i=hard-over) 

jamPos = 60; %hard-over actuator position (unit) 

setActuators (CileName, ' 05 ' , [delay;k;w;daap;dHaxHi;dHaxLo;ddotHax;hardOver; jamPos] ) ; 


^FAILED Right All Moving Tip (2nd Order, with limits); G(sj = k T w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; *;gain 

w = 63.2455532; ^freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dHaxHi = 60; *max position, upper limit 

dMaxLo = 0; ^max position, lower limit 

ddotHax = 150; %max rate (+/- unit/sec), enter positive number 

hardOver ■ 0; %hard-over switch (0=no hard-over, l=hard-over) 

jamPos = 0; %hard-over actuator position (unit) 

setActuators (f ileName, ‘ 15 ' , [delay;k;u;damp;dHaxHi;dHaxLo;ddotHax;hardOver; jamPos] ) ; 

^FAILED Left Spoiler (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; *gain 

w = 63.2455532; ^freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dHaxHi ■ 60; %ma>: position, upper limit 

dHaxLo = 0; <max position, lower limit 

ddotHax = 10; *max rate (+/- unit/3ec), enter positive number (Norn: 150) 

hardOver = 0; %hard-over switch (0=no hard-over, l=hard-over ) 

jamPos = 0; %hard-over actuator position (unit) 

setActuators (f ileName, 1 09 1 , [delay;k;u; damp ;dMaxHi ; dHaxLo ;ddotHax;hard0ver; jamPos] ) ; 
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^FAILED Right Spoiler (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s +• w A 2) 
delay * 0.00001; ^Transport delay (sec) 

k - 1; *gain 

w = 63.2455532; *freqency (rad/s) 

damp = 1.106797181; ^damping ratio 

dMaxHi = 60; *max position, upper limit 

dMaxLo - 0; km ax position, lower limit 

ddotMax » 150; imax rate (+/- mit/aec), enter positive number 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos = 0; ^-hard-over actuator position (unit) 

setActuators (f ileName, 1 19' , [delay;k;w; damp; dMaxHi; dMaxLo; ddotMax ;hard0ver; jamPos]) ; 

^FAILED Left Leading Edge Outboard Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 4- 2*damp*w*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k = 0.5; %gain 

w - 42.22321636; kt regency (rad/s) 

damp = 1.3952987; ^damping ratio 

dMaxHi = 40; %max position, upper limit 

dMaxLo = -40; ^max position, lower limit 

ddotMax = 40; imax rate (+/- unit/sec), enter positive number 

hardOver » 1; ^hard-over switch (0=no hard-over, l=hard-over) 

jamPos ® 5; ^hard-over actuator position (unit) 

setActuators (f ileName, '02' , [delay;k;w;damp;dMaxHi;dMaxLo;ddotHax;hardOver; jamPos]) ; 


^FAILED Right Leading Edge Outboard Flap (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 


delay 

= 

0.00001; 

^Transport delay (sec) 

k 

- 

1; 

*gain 

V 

- 

42.22324636; 

\t regency (rad/3) 

damp 

= 

1.3952987; 

% damping ratio 

dMaxHi 

= 

40; 

^max position, upper limit 

dMaxLo 

* 

-40; 

^wax position, lower limit 

ddotMax 

= 

40; 

^max rate (+/- unit/sec), enter positive lumber 

hardOver 

■ 

0; 

^hard-over switch (0=no hard-over, l=hard-over) 

jamPos 

» 

0; 

^hard-over actuator position (unit) 


setActuators (fileName, ' 12 ’ , [delay;k;w;damp;dHaxHi ;dHaxLo;ddotHax;hardOver; jamPos]) ; 

^FAILED Pitch Nozzle {2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*v*s + w A 2) 
delay = 0.00001; ^Transport delay (sec) 

k =1; Tfgain 

w - 39.18905204/3.918905204; J ;£reqency (rad/a) 

damp = 1.001376608; ^damping ratio 

dHaxHi = 15; ^max position, upper limit 

dMaxLo - -15; ^aax position, lower limit 

ddotMax » 10; ^max rate (+/- imit/sec), enter positive number (Worn: 60) 

hardOver = 0; ^hard-over switch (0=no hard-over, l=hard-over) NOTE; INTEGER 

jamPos * 0; ^hard-over actuator position (unit) 

setActuators (f ileName, f 10 1 , [del®y;k;w;damp;dMaxHi;dMaxLo;ddotHox;hardOver; jamPos]) ; 
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^FAILED Yaw Nozzle (2nd Order, with limits); G(s) = k * w A 2/(s A 2 + 2*damp*w*s + w A 2) 


delay 

k 

u 

damp 

dHaxHi 

dHaxLo 

ddotHax 

hardOver 

jamPos 


0.00001; ^Transport delay (sec) 

1; ^gain 

39.10905204/3.918905204; %freqency (rad/s) 

1.001376608; % damping ratio 
15; Vmax position, upper limit 

-15; <ma>: position, lower limit 

60; *max rate (+/- unit/sec), enter positive number 

Q. <hard-over switch (0=no hard-over, l=hard-over) NOTE: INTEGER 

0; ^hard-over actuator position (unit) 


setActuato t s ( f ileName , ‘ 20 1 , [ delay ;k ;w ; damp; dHaxHi ; dHaxLo ; ddotHax ; hardOver , j amPos ] ) , 


set_param(strcat (fileName , '/Actuators' ) , 'BackgrcundColor * , * red' ) 

^Initialize the FAILED state space model and controil allocation matrix 

get_param (street { fileName , ' /ICEstdl ' ) , ' A ' ,mat2str (sysFail . a,6) , 1 B ' ,»at2str ( sysFail .b,6) , C , 
mat2str (sysFail . c, 6) , 'D ' ,mat2str (sysFail . d, 6) ) ; 

^in the above line do not 3et Initial conditions becasue we want to have continuous states during the simulation 
setj?aram( street (fileName, 1 /ICEwdl 1 ) , 1 Back groundColor ' , 1 red ' ) ; 


sound(l: 100) ; *mafce a 1 second lOOhz 3ound to indicate the model was failed.... 
disp ( 1 IceFunc ti ons» Model Failed successfully!') 

errFiag = 0; 


-t 

h Do Not Fail Do Not Fail Do Not Fail Do Not Fail 

k 

case 3 

errFiag = 0; 

k 

* Case 4 Case 4 Ca3e 4 Case 4 Case 4 Case 4 Case 4 

k 

^user definable function:... call as IceFunctionsDC ( [ 1 xl ,x2 x3. . . ] ) ; 
case 4 

errFiag = 0; 


k Case 5 Case 5 Ca3e 5 Case 5 

k 

case 5 

errFiag - 0; 


Case 5 Case 5 Case 5 Case 5 


^ Case 6 Case 6 Case 6 

k 

case 6 

errFiag = 0; 


Case 6 Case 6 Ca3e 6 Case 6 Case 6 


k End Select End Select 


End Select End Select End Select End Select 


end TtEnd of select case 

^disp (' Finished IceFunctions ' ) 
return; ^IceFunctiionDC.m 
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JjsiB BBaaB CB3aE=2SSBBBDOB3 *S3n3«tSCBBCSaan = = r= CUSESSS ESSSSSCStt SB = == = == aBESa = aa = C = = lt«!E2 = = I 

function setActuators (fileName, actuatorName, values) 

^ Thia function writes the actuator properties to each actuator 
^inputs: 

Tf fileHame: from ICEfunctionsDC the name of the SIMULINK file 

h actuatorName: The name of the actuator as defined in iceDC 

^ values: The values neccessary to define the actuator 

* outputs: 

Tt Writes to a actuator in filename given. Block containing actuators must be named 

\ 'Actuators' 

^Set the parameters for a second-order acuatot model in SIMULINK 

set_j>aram (streat (fileName, ' /Actuator3/Delay 1 , actuatorName) , 'DelayTime 1 ,num2str (values (1) ) ) ; 

set_param( streat (filename, ' /Actuator 3 /Gain* ,actuatorName) , 1 Gain* ,num2str (values (2) ) ) ; 

set_param{ streat (fileName, ' /Aetuators/vw ’ , actuatorName) , 'Gain 1 ,num2str (values (3) A 2) ) ; 

set_param (streat (filename, 1 /Actuator s/d2w* , actuatorName) , 'Gain 1 ,num2str (2*values (4) *values (3) ) ) ; 

setjparaa(s treat (fileName, ' /Actuator s/d Integrator ' , actuatorName) , 'UpperSaturationLimit 1 ,num2str (values (5) ) , 

1 LowerSaturationLimit 1 ,num2str (values (6) ) ) ; 

setjmr am {streat (fileName, 1 /Actuators/ddot Integrator 1 , actuatorName) , 'UpperSaturationLimit ' ,num2str (values (7) ) , 
1 LowerSaturationLimit* ,num2str (-values (7) ) ) ; 

set_j?aram (streat (fileName, * /Actuators/Hard over Switch 1 , actuatorName) , 1 sw' ,num2str (values (8) ) ) ; 
set_paraa (streat (fileName, ' /Actuate rs/jamPositicn ' , actuatorName) , 'Value 1 ,num2str (values(9) ) ) ; 

return ^setACtuators 


function [sysFC,sysDC, sysFail,icFC, icDC]=getICEsystea (DC,FC) ; 

>ICE Linear Models 
TtUominal Model Definitions 


unemes » 

[ ’ de3 

' del3 ‘de4 

; ' deS 

' ; 'dels 

‘ ; *de9 

1 ; ' del9 

1 ; 1 ds2 

' del2 

’delO ' 

; ' de20 

’] 

xnames = 

fu 

* ; 1 w 1 ; 1 qb 

1 ; ' theta 

1 ; 1 V 

‘pb 

'rb 

' ; 1 phi 





ynames » 

[ ' ve 1 

’ ; * alphaw 1 ; ' c;b ' 

1 ; ‘ theta 

1 ; 'be taw 

'ps 

' ; *rs 

'phi 

‘axeg 

' ; ' aycg 1 

; ' aneg 

•] 


TtFirst set DC then FC 
for i = 1 : 1 : 2 , 
if i— 1; 

swtch*DC; 
elseif i«»2; 
swtch-FC; 

else 

disp( 1 error in getICEsystem ' ) ; 
end ^if 

^The following section is ldenticle to loading the selected model by choosing 'other' and typing inthe file 
%the difference is only that they are included in this file. 

^The ' t' in the matrix names indicates that it is a temporary variable which will he modified. These 
* are the exact same matrices found in the model files. 


The listing for the next section of the m-file has been omitted here, and contains the state 
space descriptions of the ICE vehicle at the various flight conditions, including the initial 
conditions on the vehicle states. The latter values are not needed in implementing the 
linear models, but are useful for describing the equilibrium flight conditions. 
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*set input command magnitudes 

the t a amp *0 ; %deg 
phiamp=5; ^deg 


otherwise *if 
k 

mdl-load(swtch) 
At « mdl.a; 

Bt = adi.b; 

Ct * mdl. c; 

Dt = mdl.d; 
uOt= mdl.uO; 
xOt* mdl.xO; 
ydt- mdl.yO; 


the text in DC or FC (depending on i) does not match any of the preloaded ca 
then read them from the file loaded and name them apropriatly. . . 


end ^select swtch 

*If you so desire to modify the ouputs ot augment any input/output matrix 
^this is the location to do so... 

At = At; 

Bt - Bt; 

Ct = Ct; 

Dt = Dt; 
uOt= uOt; 
xOt= xOt; 

xOt= zeros(8,l); VIhi 3 Sets initial conditions on the states to zero' your call 

yOt» yOt; 


* where to write the temp variables? are ue calculating DC or FC 
if i==l k then write to DC 


sysDC=ss (At,Bt / Ct,Dt ) ; 
icDC=struct( ’uO 1 ,uOt, 'xO 1 ,xOt, *y0 1 ,yOt) ; 
elseif i==2 %then write to FC 
sysFC=ss (At/Bt^Ct^Dt) ; 

icFC=struct ( 'uO l ,uOt, 'xO 1 ,xOt, 'yO' ,yOt, 'thetaamp' ,thetaamp, 'phiamp 1 ,phiamp) ; 

*Also define Failed model based on FC 
* Failed Model Definition 

TpAf ai 1 = A*l. 0Q;A£ail (4, : ) =A(4, : ) ;Afail (8, : ) «A(8, : ) ;B£ail = B ft 0. 75;C£ai 1 - C;D£ail = D; 
sysFail=sy3FC; ^initialize sysFail 
for i - 1:8 

for j = 1:8 

sysFail. a(i,} ) =3ysFC.a(i,j)* (1+0. 2* (-l) A j) ; 

, end 
end 

for i = 1:8 

for 3 = 1:11 

sysFai 1 . b ( i , j ) =sysFC .b (i, j ) * ( 1+0 . 2* ( >1 ) A j ) ; 

end 

end 

sysFail. a (4, : ) *sysFC. a(4,: ) ; ^Correct Kinematic lines back to nominal 
sysFail. a(8, : ) =sysFC. a(8,: ) ; ^Correct Kinematic line3 back to nominal 
sysFail. c = sysFC.c; 
sysFail. d * sysFC.d; 

else 

disp ( ' getICEsystem» error in the output variable definition’) 
end *if 

end k for loop 

return *getICEsystem(DC,FC) 
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E. Appendix II - Innovative Control Effector (ICE) Vehicle Description 

Introduction 

A set of linear models and trim files has been provided to aid in the development 
of a controller for an advanced tailless fighter. This document contains the variable 
definitions and units corresponding to those models. A cursory description of the 
nonlinear vehicle is offered with the decomposition of aerodynamic database. Source of 
variations in linear models, i.e., stability and control derivatives, for three different trim 
conditions is traced back to the decomposition of the aerodynamic database. Mass 
properties, vehicle geometry, eg location, and a thrust vectoring model available in the 
open literature have also been collected. The result of the linear model-derived 
reconfigurable controller will be applied to the nonlinear simulation at Langley for 


comparison with in house controllers. 

Variable List 

aerods 

dimensional aerodynamic force component along negative x-stability 

2 

axis, slug-ft/sec 

aerols 

dimensional aerodynamic force component along negative z- stability 
axis, slug-ft/sec 

alphaw 

angle of attack, deg 

alphdtw 

angle of attack rate, deg/sec 

altic 

altitude, ft 

amtl(DE5) 

left all moving tail position, deg 

amtr(DE15) 

right all moving tail position, deg 

anacc 

normal acceleration at accelerometer, g’s 

aneg 

normal acceleration at eg, g’s 

anegstb 

normal acceleration at eg along negative z-stability axis, g’s 

axeg 

longitudinal acceleration at eg, g’s 

ayacc 

lateral acceleration at accelerometer, g’s 

aycg 

lateral acceleration at eg, g’s 

betadtw 

sideslip angular rate, deg/sec 

betaw 

sideslip angle, deg 

edrag 

stability axis aerodynamic force coefficient along negative x axis 

clift 

stability axis aerodynamic force coefficient along negative z axis 

cpitch 

stability axis aerodynamic pitching moment coefficient 

croll 

stability axis aerodynamic rolling moment coefficient 

cside 

stability axis aerodynamic force coefficient along y axis 

cyaw 

stability axis aerodynamic yawing moment coefficient 

dirtrmx 

directional trim input 

dpnoz(DElO) 

pitch nozzle position, deg 

dynoz(DE20) 

yaw nozzle position, deg 

erthzdt 

velocity along inertial negative z-axis, ft/sec 

erthxic 

initial position along inertial x-axis, ft 

erthyic 

initial position along inertial y-axis, ft 

flapl(DE4) 

symmetric pitch flap position, deg 

gamma 

flight path angle, deg 

grsthr 

gross thrust, slug-ft/sec 2 
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lattrmx 

lateral trim input 

lefol(DE2) 

left leading-edge outboard flap position, deg 

lefor(DE12) 

right leading-edge outboard flap position, deg 

Iontrmx 

longitudinal trim input 

mach 

mach number 

netthr 

net thrust, slug-ft/sec 2 

pbdot 

body axis roll acceleration, deg/sec 2 

pbic, pb 

body axis roll rate, deg/sec 

phidot 

euler bank angular rate, deg/sec 

phiic, phi 

euler bank angle, deg 

plain 

power lever angle, deg 

ps 

stability axis roll rate, deg/sec 

psidot 

euler heading angular rate, deg/sec 

psiic, psi 

euler heading angle, deg 

qbdot 

body axis pitching acceleration, deg/sec 2 

qbic, qb 

body axis pitch rate, deg/sec 

rbdot 

body axis yaw acceleration, deg/sec 2 

rbic, rb 

body axis yaw rate, deg/sec 

rs 

stability axis yaw rate, deg/sec 

ssdl(DE9) 

left spoiler-slot deflector position, deg 

ssdr(DE19) 

right spoiler-slot deflector position, deg 

taill(DE3) 

left elevon position, deg 

tailr(DE13) 

right elevon position, deg 

thetdot 

euler pitch angular rate, deg/sec 

thetaic, theta 

euler pitch angle, deg 

udot 

acceleration along body x-axis, ft/sec 2 

uic, u 

velocity along body x-axis, ft/sec 

vcas 

calibrated airspeed, kts 

veas 

equivalent airspeed, kts 

vel (velkts) 

true airspeed, ft/sec (kts) 

vie, v 

velocity along body y-axis,ft/sec 

vtotdot 

rate of change of true airspeed, ft/sec 2 

wdot 

acceleration along body z-axis, ft/sec 2 

wic, w 

velocity along body z-axis, ft/sec 

wtic 

weight 


ICE Vehicle (101-3) 

The controller will be applied to the tailless fighter configuration developed under 
the Innovative Control Effectors (ICE) program. Configuration 101-3 of this program is 
shown in figure 1. The control effectors include elevons, pitch flaps, all moving tips, 
thrust vectoring, spoiler slot deflectors, and outboard leading edge flaps. The 
conventional control effectors are defined as the elevons, pitch flap, and leading edge 
flaps. The innovative control effectors are defined as the thrust vectoring, all moving tips 
and spoiler slot deflectors. Challenges associated with control using the all moving tips 
and spoiler slot deflectors include zero lower deflection limits, strong multi-axes effects 
and effector interactions (latter, not apparent in linear models).[l] 
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SKEWED AMT 



ElEVON 

SYMMETRIC PITCH 
FLAP 


PITCH/YAW THRUST 
VECTORING 


Figure 1. Model Layout 


The full nonlinear simulation of this vehicle remains proprietary. Lockheed, however, 
has granted NASA Langley Research Center permission to supply Prof. Hess with a set 
of linear models spanning the subsonic flight envelope for his reconfigurable control 
work. The supplemental information provided in this package concerning mass 
properties, etc., has been accumulated from the open literature [1], [2], allowing the 
researcher to track variations in stability and control derivatives to pertinent increments in 
the aerodynamic data base. The ultimate aim is to integrate the proposed control in the 
actual nonlinear simulation at NASA Langley Research Center to further validate the 
methodology and compare the results with the current in-house effort. Any activity 
towards addressing stability/performance robustness issues as well as nonlinear 
simulation implementation issues prior to NASA involvement will increase the likelihood 
of a successful implementation. 

Linear Models 

The set of linear models provided corresponds to three types of trim conditions: 
lg wings-level, steady level turn, and symmetric pull up/push down trim [3], The 
corresponding flight conditions span a range of Mach [.3-.9], altitude [15-35k ft], and 
angle of attack. These models provide an adequate set of stability and control derivative 
variations associated with subsonic flight. 

The pertinent trim variables are summarized in three MATLAB files: 
trmmapWLIG.mat, trmmap ST. mat, and trmmap PU.mat. These trim file summaries 
catalog the linear models respectively found in subdirectories lin_wng_lev_lg, 
lin_stdy_turn, and lin_pull_up. The MATLAB files within them are labeled as follows 
linwnglev] gmXhZ. mat 
l instdyJurnjnXa Y_hZ. mat 
lin _pull_up mXaYhZ.mat 

where 


X = 1 0 • mach 
Y = alphaw 
Z = altic/\000. 
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Along with the matrices of the state space representation, each file contains trim values of 
state, input, and output (xO, uO, and yO), and the variable names (xnames, unames, and 
ynames). Units are given above in the List of Variables. Actuator dynamics, position 
limits, and rate limits can be found in Table 1 . 

A few comments concerning the choice of trim conditions are in order. The 
loaded conditions of steady level turn and symmetric pull up/push down have been 
included to provide an independent range of angle of attack and Mach for a given 
altitude. This has been done to provide a better characterization of vehicle’s aerodynamic 
database. Note, in 1 g wings-level flight, for a given weight, speed, and altitude, there is 
only one angle of attack that will generate the necessary lift. 

The ability to independently vary Mach and angle of attack, however, comes with 
a cost. In lg wings-level flight, longitudinal and Iateral/directional dynamics are 
decoupled and the gravity vector orientation relative to the trimmed vehicle is fixed. In 
symmetric pull up/push down trim, the longitudinal and lateral/directional dynamics are 
still decoupled, but the relative orientation of the gravity vector is changing (nonzero 
thetdot). In the steady level turn, the relative orientation of the gravity vector is constant 
(zero phidot and thetadot) but the longitudinal and lateral/directional dynamics are 
coupled. The effect of the aerodynamic database on the respective derivatives is 
available, however, since the data supplied in this package includes the vehicle’s inertial 
properties and the corresponding constant components of angular rate and velocity. 
Dependence on Mach and angle of attack, for example, is available from both lg wings- 
level and pull up/push down trim conditions. An additional dependence on side-slip 
angle is available from the level steady turn trim conditions at equivalent Mach and angle 
of attack values. 


28 


Table 1. Actuator Dynamics 



Amplitude Limit 

Rate Limit 

Elevon (taill/DE3), (tailr/DE13): 



(40X100) 

(5 + 40X5 + 100) 

±30 deg 

1 50 deg/sec 

Symmetric Pitch Flap (pflap/DE4): 



(40X100) 

(5 + 40X^ + 100) 

±30 deg 

50 deg/sec 

Skewed All Moving Tip (amtl/DE5), (amtr/DEl 5): 


(40X100) 

(5 + 40X5 + 100) 

[0, 60] deg 

150 deg/sec 

Spoiler-Slot Deflector (ssdl/DE9), (ssdr/DE19): 


(40X100) 

(5 + 40 X 5 + IOO) 

[0, 60] deg 

150 deg/sec 

Outboard Leading Edge Flap (leflol/DE2), (leflor/DE12): 


(17. 828X100) 

(5 + 17.828X5 + 100) 

± 40 deg 

40 deg/sec 

Pitch/Yaw Thrust Vectoring (dpnoz/DElO), (dynoz/DE20): 


(37.186X41.3) 

(5 + 37.186X^ + 41.3) 

± 1 5 deg 

60 deg/sec 

Power Lever Angle (plain)* 



Afterburner: ^ 

.555 + 1 

[30, 90] deg 

22 deg/sec 

Afterburner out: ^ 

.6255 + 1 

[90, 127] deg 

14 deg/sec 


’•'dynamics currently not in proprietary simulation, but proposed for LaRC version 

Nonlinear Simulation and Linear Models 

Some features of the nonlinear simulation are considered in the context of the 
linear models generated. Components of the trimming forces and moments are discussed 
along with the impact of various aerodynamic increments on the stability and control 
derivatives for the flight conditions considered. The mass properties data and thrust- 
vectoring model are available in the open literature. Consider first a cursory look at the 
mathematical framework of the nonlinear vehicle description. 
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The ICE model is governed by the standard set of equations of motion for a flat- 
earth, rigid-body, symmetric aircraft [4], 


“ = ~% w b + r b v b -gsi”0 + ~ 

m 

v = -r b u + p b w + g cos 6 sin <j> + ¥- 

m 

w = -p b v + q b u + g cos 6 cos </> + — 

m 

Pb'b^xx - 0 + r lhz - plhz\ m 


yy 

Pb l_ 

r b 


XX ] xz 

-hz hz 


1-1 


L -% r b ( I r ~ OtPb % 1 X z\ m 

N - PbVbVyy ~ I zz)~ < ]b r b 1 xz) ^ 


180 


(la) 

(lb) 

(lc) 

(l d) 

(le) 


9 = q b cos <j>-r b sin <j> 

$ = Pb + q b tanO sintf) + r b tan9 cos<f> 


(10 

(lg) 


The model constants are listed in Table 2 [1]. The force and moment vectors, 
F =\X Y z] and M =\L M A'], respectively, expand as 


F = F a + F T 

K = K,o + K.dvn + K,S 


( 2 ) 


F T = F'+F rd 


M = M a + M t 
K = M a ' 0 +M adyn + M a , s 
M j = M t + M r( j 


( 3 ) 


where subscripts ‘ a ’ and ‘ T ’ correspond to the aerodynamic and thrust related forces 
and moments. The aerodynamic force and moments are further decomposed into ‘ a,o ’ 
components of the baseline aircraft with no effectors deflected, ‘ a.dyn ’ components due 
to angular rates, and ‘ a, 8 ’ force and moment increments due to control deflection. The 
thrust force and moments are decomposed into ‘ t ’ components due to gross thrust and 
‘ rd ’ components due to ram drag. Although it is unlikely that anything more than a 
first-order approximation of the aerodynamic data can be obtained from the trim files and 
linear models provided (maybe not with the decomposition above), it is of interest how 
these terms influence the stability and control derivatives particularly when 
characterizing uncertainty. In the coming subsections, the aerodynamic data build up is 
considered along with its influence on the linear models. 
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Table 2. Model Constants 


Symbol 

Description 

Value 

J 

jc-axis moment of inertia 

35,479 slug-ft 2 

l yy 

y-axis moment of inertia 

78,451 slug-ft 2 

yy 

i Z2 

z-axis moment of inertia 

1 10,627 slug-ft 2 

i X2 

xz-axes product of inertia 

-525 slug-ft 2 

S ref 

wing planform area 

808.6 ft 2 

c 

wing mean aerodynamic chord 

28.75 ft 

b 

wing span 

37.5 ft 

m 

Aircraft mass 

1017.9 slug 

w 

Aircraft gross weight 

32,750 lb 

“ TWL 

waterline distance between eg and thrust 
application point 

-.417 ft** 

“ TFS 

fuselage station distance between eg and 
thrust application point 

18.79 ft** 

d DWL 

waterline distance between eg and ram 
drag application point 

-0.3308 ft 

d DFS 

fuselage station distance between eg and 
ram drag application point 

-12.66 ft 

^ACWL 

waterline distance between eg and 
accelerometer 

-1.039 ft 

dpSFS 

fuselage station distance between eg and 
pilot station 

-13.55 ft 

d AC FS 

fuselage distance between eg and 
accelerometer 

-11.5008 ft 


* relative to eg position of 38% of mean aerodynamic chord, nominal weight 
** positive WL, FS distance is above and behind eg 


Aerodynamic Forces and Moments 

The total aerodynamic forces and moments at trim can be obtained from the trim 
file as follows 


X a,s=-V S ref cdra £ 

(4a) 

Y a, S = <jS refute 

(4b) 

Z a ,=-qSref ’Clift 

(4c) 

L a,s=dS ref bCTOll 

(4d) 

M a,s =qS ref ccv\\ch 

(4e) 

N a ,s = *}S ref b • cyaw . 

(4f) 
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These components include trim forces and moments due to the baseline vehicle’s 
orientation and speed with respect to the air stream, due to any trim angular velocities 
(excluding inertial effects), and due to the trim control deflections. Note these 
components are expressed with respect to stability axis. Transforming these components 
to body axis as reflected in the equations of motion is straightforward. The 
decomposition of components into the corresponding ‘ a,o \ ’ a,dyn \ ’ a, 5' terms may 

not be so straightforward or possible. Some additional background information on the 
dynamic term is provided. 

Rotary and Forced Oscillation Terms 

The direct method in AIAA-88-4357 is used to add the rotational and oscillatory 
increments to the force and moment components for a < 30° . The effect of angular rates 
on the aerodynamic moments and forces is modeled using data from two sources: 1) a 
rotary balance test where the vehicle rotates about the velocity vector and 2) a forced 
oscillation test where the vehicle is forced to oscillate about each body axis. 

First the angular velocity vector is equivalently expressed as 

® = pJb + %h + r bh = °>vK + “> os 

(o v - 6) ■ i w = p b cos a cos /? + q b sin/3 + r b sinacos J3 

where i w is a unit vector along the velocity vector, (o v and cb os are respectively the 

rotary (along the velocity vector) and oscillatory components of the angular velocity 
vector. The oscillatory component is further decomposed along the body axes as 

“>os = PoJb + %s)b + r osK ( 6 ) 


where 


Pos = Pb~ cos a cos P (7a) 

=q b -co v sinp (7b) 

r os = r b~ My sin a cos/}. (7c) 

From these expressions, it is clear that this decomposition provides an equivalent 
representation of the angular velocity. The components of the oscillatory term and the 
rotary angular velocity are made non-dimensional as follows 


CO os = 


hPo 

2V 


II + 


c %i 

2V 


Jb + 


b _L k 

2V 


(Oy = 


bcoy 

2V 


(8a) 

(8b) 


The direct method uses the non-dimensional terms to define dynamic effect of angular 
rates on the forces and moments. 


P a,dyn F aro {M ,CC, f 3 ,G)y} ( 9 a) 

M a ,jy n = M a jM,a,/1^)+njM,a)a>’ os (%) 


32 



Subscripts ‘ a, ro ’ and ‘ os ’ denote respectively the rotary and oscillatory components. 
The matrix FI M consists of the following oscillatory dynamic derivatives 


n„ 


Lp,os ® L r os 

0 0 
N p „ 0 N'„ 


( 10 ) 


The rotary term is zero whenever co v is zero, a result of either 6) = 0 or d> is 
perpendicular to I w . The former occurs in lg wing-level trim and the latter occurs in 

both the symmetric pull up/push down and steady level turn. Consequently, the trim 
values of force and moments do not have any rotary component. The rotary components 
do influence the derivatives. It should be mentioned that the Mach dependence in the 
rotary terms exclusively resides with dynamic pressure, q , used to dimensional-ize the 
corresponding force and moment coefficients. The axial force rotary component is zero. 

At trim the steady effects due to constant angular rates, if any, result from the 
oscillatory component, or 


F =0 

r a ,dyn, trim u 

^o,dyn, t nm=^os{ M ’ a Ks- 


(11a) 

(lib) 


The aerodynamic trim force and moment components, then, take on the following form 
1- g, wings level: p os = q’ os = r' os = p = ° 


K 1 r, m = Ko( M ' a 'fi) + KA M - a ’ s ) 

M ur , m =M a jM.a,p) + M aS ( M.a.S ) 
pull up/push dow n: p'os = r' os = p = ° , q' os = cq b j 2V 


a,trim 


K.o(M,a,p)+F aS (M,a,S) 




0 

M q.os^os 

0 


+ M atS {M,a,S ) 


(12a) 

(12b) 


(13a) 

(13b) 


steady level turn: p' os = bp b /2V , q' os = cq b j 2V , r' os = br b /2V , /? = constant 


K, m =P>, a ,p)+KA M ' a ' s ) 


K Jrlm =K :0 (M.a,Ph 


Lp'OsPos ^ A-.as^as 
^q.osQos 

^ p.osPos os^OS 


+ M a6 (M,a,S). 


(14a) 

(14b) 


If a first-order approximation is assumed for F aS and M a S about S = 0 (trim effector 
deflections are small), so 
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(15a) 

(15b) 


Ks = Kjf* 

M a , s = K,s d . 


F ao and M ao can be approximated for lg-level flight. Note, F aS and M .' aS are 
available from the B matrix. Altitude effects can be removed by normalizing with q . It 
will be shown that for both lg-level flight and symmetric pull up, the term M qos is 

explicit in the A-matrix. As a result, F a 0 and M ao can be approximated for symmetric 
pull up. Side-slip angle effects, only on F ao , can be determined from the difference 
obtained in the term for the symmetric pull up and steady level turn at the same (a,M ) 
flight conditions. Side-slip angle effects on M ao (except for the pitching component) as 

well as identifying the remaining dynamic derivatives may not, however, be possible. 
Stability Derivatives and the Nonlinear Aerodynamic Database 

Let us assume that A' is the linear system’s A-matrix with all the differential 

contributions due to inertial terms (e.g., wq b , uqhJ^p^, [l zz ~ lyy^b^ ^ gravity 

terms removed. Only the aerodynamic contribution remains in A'. Pursuing a 
construction similar to that use in [4], the differential forces and moments determine the 
linear models, i.e., what is leftover in A' 


a 88 


dM„ 


88 


d -^dc 

da 

jK 

da 


dF 


dF 


dF 

dp b + 

dF, 

J dF a , 

+ — - 

dp 

+ — - 

dM- 

a 


dq b + ® dr b 

d(i 

dM 


d Pb 


d % 

dr b 

dM 


dM 

dM h 

dM n 

d P b 

dM . 

da + - 

a 

dp + 

a 

b a 

+ , d d b + 


dp 

dM 


d Pb 

Sq b 


dr. 


where all partial derivatives are evaluated at the trim condition. The terms on the right 
influence the A-matrix whereas the second term (negative of) on the left corresponds to 
the B-matrix. Substituting the differentials da , dp , and dM for expressions in table 3 
expressed in terms of du , dv , and dw provide the derivatives used to construct A. Note 
when the trim /? (v) is zero as in 1-g level flight or symmetric pull up, dM influences 
only the du dw columns of A, whereas dp influences only the dv column of A. In a 
steady level turn, dM and dp influence all columns of A corresponding to translational 
velocity. From its definition, da influences only the du dw columns of A regardless of 
the trim. 

It should be mentioned that the columns corresponding to du and dw are also 
influenced by the differential force and moments due to thrust, dFj and dMj , since the 
application of ram drag is dependent on angle of attack as will be shown later. 
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Table 3. Transformation to translational velocity 



dw 


*all values at trim; a is the speed of sound. 

The contribution of the baseline aerodynamics and the increment due to control 
deflection 




AK.C ±Ks) da | d im d p + F'„„ + K;> ) 


da 


dp 


dM 


dM 


i - - \ dM nX 


dfe,., + ^<o d/1 + #«L + Mod 


da 


dp 


dM 


dM 

(17) 


does not influence the columns of A corresponding to dp b , dq b , and dr b . The dynamic 
aerodynamic increment 


r)F r)F 

dF a .=-^da + -^- 
a4yn da dp 

dK'dyn SM 


dp + 


dF, 


a,dyn 


dM 


a t dyn 


da 


da + 


a,dyn 


dM 
dM 


dM + - 


dF, 


a,dyn 


d Pb 


d Pb+- 


dF„ 


a t dyn 


dp 


dp + a J^dM + 

dM dp b 


d % 

d Pb + 


d % + 


dF r 


a,dyn 

drF 


dr. 


dM 


a,dyn 


dq h 


d % + 


dM 


a,dyn 


dr, 


dr. 


b 

(18) 


potentially affects all six columns. Next, expansions of the partials in equation 18 are 
developed. The result will be the dynamic aerodynamic increment has no influence on 
the derivatives corresponding to translational velocities in lg-level flight. The influence 
of the dynamic term expands, however, as trim pitch rate becomes nonzero for symmetric 
pull up flight, and the rest of angular rate components and side-slip angle become 
nonzero in the level steady turn. 

For partial derivatives taken with respect to jc ( = a,p, or M , the initial form is the 

same 


dF, 


a,dyn 


OX: 


dF a,ro + dF a,ro do'y 


dX; 


d(o' v dx ; 


dM a, dyn _ + dM aro dcpy + 


dx. 


dx, 


dco\, dx ; 


— ^ cb ' 
dx- os 


+a 


, dd) 'os day 


dXj da>y dx i 


(19a) 

(19b) 
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Since the rotary component is zero whenever co' v is zero (true for all trim), the first term 
in (19a-b) is always zero. The second term is zero for = M in both equations, since 

M = ( 20) 

dM V v 


The third term in (19b) is zero when x t = . The components of (19a-b) can be found in 

tables 4 and 5, where the following shorthand is used. L l0 <^ ro , M a ^ 0 , X M ^ ro are the 

components of dM aro /do^; X <ro ,Y <ro ,Z <ro are the components of dF aro /d(o' v ; 
Lposa anci LposM are the partials of L p os with respect to a and M , respectively; s a 
and Cp are sin(a) and cos(j3), respectively; stability axis roll and yaw rate are denoted 
as p s and r s , respectively. 

Table 4. Dynamic Force Partial Derivatives with respect to (a, J3,M) 


^_dyn _ q 

da 


^ Y dyn _ _b_y L c ) 

da 2V P ’ 


dZ dy n = J>_z ( rc ) 

da 2V P ’ 



dp 

Mjfyn 

dp 

dZ dyn 

dP 


^ Z ^ C P- P S S P) 


dX. 


dyn 


dM 


dY, 


dyn 


dM 


dZ. 


dyn 


dM 


= 0 
= 0 
= 0 
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Table 5. Dynamic Moment Partial Derivatives with respect to (a, /?, M) 


5 

da 

dM 


a 2V ^«> v ,ro sC P + L P ° s a Pb + Lr -°s a r b r s(j'p,os C a C ft + K,os S a C ft Jj 


— = — M , r s c g +-£-\ M a osa q b - r -j-M a os s 20 

da 2V (0 v r0 P 2V l Q’ os ' a ^° 2 q ' 


dN 


= Tf. 

da 2V V c °v ro P 


+ N p.os a Pb + N r,os a r b - r s{ N p,os C a C ft + N r,os S a C ft] 


dp 

dM 


'j- = w[ L <ro (WP ~ P s S pV{ L P,os{\ Ps C a S 2ft ~ % C c/ft) + L r,os{~ Ps S a S 2ft ~ Wft 


dp 


— M co v ,ro (wp -Ps S ftV~{ M q .o[p/p ~ \ %hp 


dN 


dyn _ _b_ 


dp 2V 


{ N <r jpb C ft ~ Ps S ft) + PM P ~ ^a^p] + N r, ( ,( 1 2 Ps S a S 2ft ~ W 


dL 




dM 2V 
dM 


■' p,os,M 


y Lp.os^Pb i^r,os,M y 


dM 2V\ q,osM V Mq ' os 


dN dyn _ b ( ( 

dM 2V 


N P.os,M ~y N p,os \Pb +| N r,os.M ~~ N r.os fb 


In 1 -g level flight, all terms in tables 4 & 5 vanish due a zero angular velocity vector at 
trim. For symmetric pull up, there are nonzero entries for dM adyn jdp and the 
components dM dyn jda and dM dyn jd M due to the nonzero trim q h . For level steady 

turn, all expressions not specified as zero are nonzero. 

For partial derivatives of equation 18 taken with respect to angular velocity, the 

form is 


d Kdyn _ dP a,ro d(0 V 

deb dco'y deb 


a,dyn _ dp)y ^ j-j 


d(b 


deo'y 


deb 


r <Ks , d “>'os d(0 V X 


deb 


deoy 


deb 


(21a) 

(21b) 


Tables 6 and 7 contain the nine components associated with each equation. The 
shorthand used in tables 4 and 5 is again used here. 
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To be noted, when /? = 0 as in 1 g-wings level and symmetric pull up flight, all partials 
with respect to q b vanish except 


<M_ = m 

dq b 2V 


q,os ' 


( 22 ) 


38 


For these two cases, the longitudinal pitch rate (perturbation) does not effect the 
lateral/directional moments. However, the lateral/directional roll and yaw rate 
(perturbations) do influence (albeit small from observation) the pitching moment through 
the rotary term. It appears that both the rotary and oscillatory components influence the 
lateral/directional moments in response to perturbations of the lateral/directional angular 
rate components. In the steady level turn, the longitudinal pitch rate (perturbation) does 
effect the lateral/directional moments providing some of the aerodynamic coupling at this 
trim condition. The rest of the aerodynamic coupling arises from the expressions in 
tables 3 and 4. 

To close this discussion on the source of the linear models, the thrust-vectoring 
model is presented followed by the trim control used. 

Thrust Equations 

Thrust vectoring is applied to both the pitch and yaw axes. The thrust force and 
moment components may be expressed as: 


X t = T g COS { 5 <p) COS { 5 ty) 

(23a) 

Y t =T g cos\$ tp )sin{8 t y) 

(23b) 

Z t =-T g sin{5 tp ) 

(23c) 

II 

1 

t'-T'"* 

(23d) 

M t =l x Z l+ l z X t 

(23e) 

= -IJ, 

(230 

where T g is gross thrust (lb), 8 tp (dpnoz) and 5 ty 

(dynoz) are the pitch and yaw nozzle 

deflections, respectively, and where l z and / are 

: respectively the distances (ft) below 

and behind the eg for the thrust application point. 
T g =grsjhr; l z = -d Tm ; l x = d TFS . 

In terms of variables defined above. 

Force and moments due to ram drag, D ram (lb), are expressed (best guess) as 

1 

II 

Q 

ram g nei 

(24) 

X rd = ~ D ram C0S ( a ) 

(25a) 

o 

II 

(25b) 

Z rd ~ ~ D ram si 4 a ) 

(25c) 

o 

II 

(25d) 

^ rd ^x.rd^rd ^z,rd^ rd 

(25e) 

O 

II 

s? 

(250 


where T m( is net thrust, and where l x rd and l z rd are the distances(ft) ahead and below 
eg for the ram drag application point. From the trim file and Table 1 .0, T net = net thr; 

( x.rd ~ ~^DFS ’ h.rd ~ ~^DWL • 
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The differential force and moments due to thrust are used in linear models. 


H 

dS tn tp dd, 




ty 


dF T = dF t + dF rd = 


+ u / \ Jf 

dM r = dM, + dM, = — ± dS ln + -— *■ dS w + 

T 1 rd 8S tp tp d5, y ty da 


iy 


da 


dM, 


r(L d< 


( 26 ) 

(27) 


Equations 23 and 25 determine the partials above. The first two partial derivatives of 
equations 26 and 27 determine the columns of B corresponding to thrust vectoring. The 
third modifies the du,dw columns of the A-matrix. 

Trim control used in linear models 

Lastly, three inputs, ‘lontrim’, ‘lattrim’, and ‘dirtrim’ drive the control effectors to 
respectively trim the vehicle about the longitudinal, lateral, and directional axes. 

elevl - lontrim + lattrim 
elevr - lontrim - lattrim 
pflapx = 0. 1 1 1 1 • lontrim 

leflox = 0. 07 1 1 1 • [lontrim + lattrim) 
lefrox = 0. 07 1 1 1 (lontrim - lattrim) 
dpnozx = 0. 1 6 • lontrim 
dynozx = 0. 1 6 • dirtrim 
ssdl = ssdr = amtl = amtr = 0 


For each effector, the weighting used is the squared value of the effector’s rate limit 
normalized with the elevon’s rate limit. Due to difficulties with the trimming routine, the 
unilateral controls (those restricted to positive deflections) were not used in the trim. The 
linear models for the four pertinent flight conditions are presented next. 


Linear Models 

State space representations for linearized ICE models for the flight conditions of interest 
for the design challenge are given on the following pages. The state, output and input 
vectors in these linear models are defined as: 


x T = [u w qb theta v pb rb phi] 
y = [vel alphaw qb theta betaw ps rs phi axcg aycg azcg] 
u=[dE3 dE13 dE4 dE5 dE15 dE9 dE19 dE2 dE12dE10 dE20] 
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Mach No. = 0.3, Altitude = 15,000 ft (wings level) 


5.9570e-003 
-6. 3377e-O02 
-5. 2181e-002 
0 

3. 32X6e-007 
8. 5005e-006 
-2. 9015e-006 
0 


3 . 1987e-002 
-5. 8625e-00i 
8. 1704e-002 
0 

8 . 2150e-008 
2 . 1012e-006 
-7 . 1571e-007 
0 


-1 . 3310e+000 
5. 3737e+000 
-5 . 9437e-001 
1 . OOO0e+OO0 
0 
0 
0 
0 


-5. 4507e- 
-1 . 3501e- 


001 

001 

0 

0 

0 

0 

0 

0 


-1 . 2579e-003 000 

-1.9032e-003 2.3209e-004 5.7485e-005 0 

1 . 1259e-002 1.7539e-003 4.3442e-004 0 

0 0 0 0 

1. 1094e-002 1.3318e+000 -5.3735e+000 5.4507e-001 

-1 . 2207e+000 -6.0701e-001 3.9515e-001 0 

-1. 2626e-Q01 2.9703e-003 2.0646e-002 0 

0 1. 0000e+000 2 . 4768e-001 0 


» b 


b = 


Columns 1 through 8 

-5. 9693e-003 -5.9693e-003 -4.0856e-003 
-2. 6797e-001 -2.6797e-001 -2.3059e-001 
-1. 3850e+000 -1.3850e+000 -1.2023e+000 
ooo 
-9. 1236e-003 9.1236e-003 0 

2.4117e+000 -2.4117e+000 0 

3. 7657e-002 -3.7657e-002 0 

0 0 0 


-1 . 4174e-002 -1.4174e-002 
-8 . 4039e-002 -8.4039e-002 
-3. 5207e-001 -3. 5207e-*001 


-2 . 4849e-004 
5. 3690e-001 
-5. 6093e-002 
0 


-2. 5238e-002 
1. 1122e-001 
2. 6650e**001 
0 0 

2. 4849e-004 2.1933e-002 

5. 3690e-001 -1.1963e+000 
5.6092e-002 -1.2350e-001 
0 0 


-2. 5238e-002 
1.1122e-001 
2. 6650e-001 
0 

-2. 1933e-002 
1. 1963e+000 
1 . 2350e-001 
0 


1 . 1500e-002 
9. 7043e-002 
2. 3341e-001 
0 

-2 . 8901e-002 
-9. 8532e-001 
1 . 8590e-001 
0 


Columns 9 through 11 


1. 1500e-002 -1 . 1805e-005 
9. 7043e-002 -1.2395e-001 


2. 3341e-001 
0 

2. 8901e-002 
9. 8532e-001 
-1. 8590e-001 
0 


-1. 7312e+000 
0 
0 
0 
0 
0 


0 

0 

0 

0 

1 . 2395e-001 
-6. 6733e-002 
-1 . 2273e+000 
0 


» c 


9. 7067e-001 
-4. 3428e-002 
0 
0 
0 
0 
0 
0 

1 . 8515e-004 
1. 0324e-008 
1. 9698e-003 


2. 4042e-001 
1. 7534e-001 
0 
0 
0 
0 
0 
0 

9 . 9420e-004 
2. 5533e-009 
1. 8221e-002 


0 

0 

1. 0000e+000 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

1 . 0000e+000 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

1.8063e-001 

0 

0 

0 

-3. 9096e-005 
3. 4482e-004 
5. 91S3e-005 


0 

0 

0 

0 

0 

9. 7067e-001 
■2. 4042e-001 
0 
0 

2. 7166e-005 
-7. 2136e-006 


0 

0 

0 

0 

0 

2. 4042e-001 

9. 7067e-001 
0 
0 

6. 7286e-006 
-1 . 7867e-006 


0 

0 

0 

0 

0 

0 

0 

1 . 0000e+000 
0 
0 
0 
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» d 


Columns 1 through 8 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1. 8553e-004 

-1. 8553e-004 

-1 . 2698e-004 

-4.4053e-004 

2.8357C-004 

2. 8357e-004 

0 

-7. 7233e-006 

8. 3287e-003 

8. 3287e-003 

7.1669&-003 

2. 6120e-003 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-4. 4053e-004 -7.8444e-004 

-7. 8444e-004 

3. 5743e-004 

7. 7233e-006 

6.8170e-004 

-6. 8170e-004 

-8. 9826e-004 

2. 6120e-003 -3.4567e-003 

-3. 4567e-003 

-3. 0162e-003 


Columns 9 through 11 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

ooo 
ooo 
ooo 
ooo 
3. 5743e-004 -3.6690e-007 0 

8. 9826e-004 0 3.8524e-003 

-3. 0162e-003 3.8524e-003 0 
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Mach No. = 0.6, Altitude = 25,000 ft (wings level) 


-4.4568e-003 
-5.2647e-002 
-1. 0290e-002 


6.4410e-002 
-5 . 5615e-001 
-1 . 7800e-001 


0 

0 

1. 5401e-019 
-1 . 0408e-017 
0 


0 

4. 2110e-018 
4. 6198e-019 
-3.1220e-017 
0 


-1. 1655e+Q00 
1. 0575e+001 
-5. 8883e-001 
1.0000e+000 
0 
0 
0 
0 


-5. 5816e-001 -3. 2481e-003 
-6 . 1516e-002 1.2911e-002 


2. 8958e-002 
0 

-3. 4652e-003 
-8. 8575e-001 
-1 . 4300e-001 
0 


0 

-9. 1850e-005 
-3. 7347e-003 
0 

1. 1675e+000 
-8. 0642e-001 
7. 6043e-003 
1. 0000e+000 


0 

-1. 0123e-005 
-4. H61e-004 
0 

-1. Q575e+001 
1. 7919e-001 
-6. 1956e-003 
1. 1021e-001 


0 

0 

0 

0 

5 . 5816e-001 
0 
0 
0 


» b 


Columns 1 through 8 


-9. 4029e-003 
-5. 9531e-001 
-4. 2251e+000 
0 

4. 5640e-003 
5. 9008e+000 
-4. 2362e-002 
0 


-9. 4029e-003 
-5.9531e-00i 
-4. 2251e+000 
0 

-4. 5640e-003 
-5 . 9008e+000 
4. 2362e-002 
0 


-2. 7374e-002 
-3. 0032e-001 
-2. 048ie+000 
0 
0 
0 
0 
0 


-1 . 7097e-002 -1. 7097e-002 


-3. 4040e-001 -3.4040e-001 
-1 . 6221&+000 -1 . 6221e+000 
0 0 

9 . 6242e-003 -9.6242e-003 
2 . 5672e+000 -2 * 5673e+000 
-1 . 5954e-001 1.5958e-001 

0 0 


-9. 7633e-002 
2. 3746e-001 
7. 9212e-001 
0 

5. 6742e-002 
-2. 7252e+000 
-4. 7640e-001 
0 


-9. 7633e-002 
2. 3746e-001 
7. 9212e-001 
0 

-5. 6742e-002 
2. 7252e+000 
4. 7640e-001 
0 


-8. 3672e-003 
4. 6687e-003 
-1 . 0211e-001 
0 

1 . 2412e-002 
-6. 1473e-002 
-6. 8770e-002 
0 


Columns 9 through 11 


-8. 3672e-003 -8.4481e-005 
4. 6687e-003 -7.7582e-002 
-1 . 0211e-001 -1 . 0836e+000 


0 0 

-1 . 2412e-002 0 

6. 1473e-002 0 

6. 8770e-002 0 

0 0 


0 

0 

0 

0 

7. 7582e-002 
-4. 1771e-002 
-7. 6822e-001 
0 
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» c 


9. 9398e-001 
-1 . 0297e-002 
0 
0 
0 
0 
0 
0 

-1.3852e-004 

0 

1 . 6363e-003 


I. 0955e-001 
9. 3426e-002 
0 
0 
0 
0 
0 
0 

2 . 0019e-003 
1 . 3088e-019 
1 . 7286e-002 


0 

0 

1 . OOOOe+OOO 
0 
0 
0 
0 
0 
0 
0 
0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1. OOOOe+OOO 

0 

0 

0 

9. 3992e-002 

0 

0 

0 

9. 9398e-001 

0 

0 

-1. 0955e-001 

0 

0 

0 

0 

-1.0095e-004 

0 

0 

-1. 0770e-004 

6. 2015e-005 

0 

-4. 0128e-004 

2. 8548e-006 


0 

0 

0 

0 

0 

1 . 0955e-001 
9. 9398e-001 
0 
0 

6. 8348e-006 
3. 1463e-007 


0 

0 

0 

0 

0 

0 

0 

1 . OOOOe+OOO 
0 
0 
0 


» d 
d * 


Columns 1 through 8 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


-2. 9225e-004 -2. 9225e-004 
1 . 4185e-004 -1.4185e-004 
1 . 8503e-002 1.8503e-002 


0 

0 

0 

0 

0 

0 

0 

0 

-8. 5080e-004 
0 

9. 3342e-003 


-5. 3140e-004 -5. 3140e-004 
2. 9913e-004 -2. 9913e-004 
1. 0580e-002 1.0580e-002 


0 

0 

0 

0 

0 

0 

0 

0 

-3.0345e-003 
1 . 7636e-003 
-7. 3805e-003 


0 

0 

0 

0 

0 

0 

0 

0 

-3. 0345e-003 
-1. 7636e-003 
-7. 3805e-003 


0 

0 

0 

0 

0 

0 

0 

0 

-2. 6006e-004 
3. 8579e-004 
-1. 4511e-004 


Columns 9 through 11 


0 

0 

0 

0 

0 

0 

0 

0 

-2. 6006e-004 
-3. 8579e-004 
-1 . 4511e-004 


0 

0 

0 

0 

0 

0 

0 

0 

-2. 6258e-006 
0 

2. 4113e-003 


0 

0 

0 

0 

0 

0 

0 

0 

0 

2. 4113e-003 
0 
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Mach No. = 0.6, Altitude = 25,000 ft (steady turn, bank angle = 60 deg) 


8.8389C-003 

1.4825e-002 

8.4736e-002 

0 

-4.4602C-002 
-2. 1457e-003 
-6.9972e-004 
0 


-4. 5946e-002 -2.0301e+000 -5.5886e-001 
-8.5038e-001 1 . 0444e+001 -2.7508e-002 

-7. 5301e-001 -7.2466e-001 0 

0 5.0158e-00X 0 

-8. 6740e-003 1.1329e-D06 -4.7444e-002 

6. 4421e-003 -4. 0192e-002 0 

2. 5245e-003 3.8442e-D03 0 

0 8 . 4895e-002 8.9769e-002 


3. 9864e-002 0 

-8. 9107e-003 -8.0981e-003 
-1.3201e-001 3.6525e-002 

0 0 

1.1079C-002 2.0316e+000 

-1.4822e+000 -8. 8882e-001 
-2.3373e-001 -2. 3010e-002 
0 1. OOOOe+OOO 


» b 


Coluans 1 through 10 


3. 9019e-002 3.9019e-002 

-5.0400C-001 -5. 0400e-001 
-3. 9298e+000 -3.9298e+000 
0 0 

4.5591e-003 -A. 5591e-003 
5. 7041e+000 -5.7041e+000 
1. 0613e-001 -1 . 0613e-001 
0 0 


3. 7961e-002 -2.7021e-002 -2.7039e-002 -8.0999e-002 
-3. 0129e-001 -2.0913e-001 -2.0889e-001 2.3251e-001 

-1 . 9451e+00D -1.0485e+000 -1.0468e+000 6.9107e-001 

0 0 0 0 
0 1. 4190e-002 -1.4249e-002 5.6379e-002 

0 1. 5498e+000 -1. 5444e+D00 -3.0085e+000 

0 -2. 1517e-001 2.1538e-001 -3.9858e-001 

0 0 0 0 


Column 11 

7. 6809e-005 
0 

2. 3792e-005 
0 

1.9504e-001 

-1.0501e-001 

-1.9313e+000 

0 


7. 6618e-003 
-8. 4805e-005 
-1. 0136e-002 
-8.6511e-001 
-1. 0443e+001 
2. 9510e-001 
2. 2584e-002 
4. 9222e-002 


0 

-4. 8347e-001 
0 

-8. 8912e-002 
2. 8031e-001 
0 
0 

-8. 7771e-016 


-8. 1138e-002 
2. 3314e-001 
6. 9652e-001 
0 

-5. 6504e-002 
3.0183e+000 
3. 9945e-001 
0 


-5. 3D58e-003 
2.2552e-002 
-9. 9962e-002 
0 

2. 3896e-003 
-3. 1068e-001 
-1. 3590e-002 
0 


-5 . 3058e-003 
2.2552C-002 
-9. 9962e-002 
0 

-2.3896e-003 
3. 1068c-001 
1 . 359Oe-0O2 
0 


9. 3231e-004 
-1.9504e-001 
-2.7238e+0D0 
0 

-3. 6715e-007 
1.976Be-007 
3 . 6355e-006 
0 


» C 


9 . 8163e-001 

1. 9081e-001 

0 

1 . 7934e-002 

9. 2265e-002 

0 

0 

0 

1. OOOOe+OOO 

0 

0 

0 

6. 6444e-005 

-1. 2915e-005 

0 

8. 1498e-004 

4. 1927e-003 

0 

9. 9344e-007 

5. 1108e-006 

0 

0 

0 

0 

2. 7472e-004 

9. 6267e-004 

0 

1 . 5245e-007 

1 . 5897e-0D6 

3. 5211e-008 

1 . 9299e-003 

2 . 6431e-002 

9. 9481e-009 


0 

7. 2015e-004 

0 

0 

0 

0 

0 

0 

0 

1. OOOOe+OOO 

0 

0 

0 

9.3992e-002 

0 

0 

0 

9. 8163e-001 

0 

0 

-1. 9081e-001 

0 

0 

0 

0 

-1. 4710e-004 

0 

0 

3. 4434e-004 

4.7996e-005 

0 

5. 4814e-004 

1. 3560e-005 


0 

0 

0 

0 

0 

1. 9081e-001 
9. 8163e-001 
0 
0 

9. 3295e-006 
2. 6358e-006 


0 

0 

0 

0 

0 

0 

□ 

1. OOOOe+OOO 
0 
0 
0 
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>> a 


Columns 1 through 10 


0 

0 

0 

0 

0 

0 

0 

0 

1 . 2128e-003 
1 . 4170e-004 
1.5665e-002 


0 

0 

0 

0 

0 

0 

0 

0 

1. 2128e-003 
-1. 4170e-004 
1 . 5665e-002 


0 

0 

0 

0 

0 

0 

0 

0 

1.1799e-003 

0 

9.3644e-003 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
ooo 
0 0 0 
0 0 0 
0 0 0 
-8. 3983e-004 -8.4040e-004 -2.5175e-003 
4.4105e-004 -4. 4286e-004 1.7523e-003 

6 . 5000c-003 6. 4924«-003 -7.2266e-003 


0 

0 

0 

0 

0 

0 

0 

0 

-2. 5218e-003 
-1 . 7562e-003 
-7. 2461e-003 


ooo 

ooo 

ooo 

ooo 

ooo 

ooo 

ooo 

0 D 0 

-1.6491C-004 -1 . 6491e-004 2.8977e-005 

7.4271C-005 -7.4271C-005 -1.1411e-008 
-7.0093C-004 -7. 0093e-004 6.0621e-003 


Column 11 


0 

0 

0 

0 

0 

0 

0 


2. 3873e-006 
6.0621C-003 
0 
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Mach No. = 0.9 Altitude = 35,000 ft (wings level) 


1. 1400e-002 
7.1202e-002 
2.6984e-001 
0 

1. 2015e-017 
0 
0 
0 


4. 9835e-002 
-6. 1261e-001 
-4. 5007e-001 
0 
0 
0 
0 
0 


-1 . 1490e+000 -5. 5995e-001 
1. 5239e+001 -4. 2221e-002 


-6. 0733e-001 0 

1 . 0000e+000 0 

0 0 

0 0 

0 0 

0 0 


-4. 3042e-003 
8 . 9952e-003 
2.4971e-002 
0 

-8. 1296e-003 
-7. S942e-0D1 
-1 . 2068e-001 
0 


0 

3. 9819e-004 
-7. 0499e-004 
0 

1. 1506e+000 
-7. 5934e-001 
8, 4167e-003 
1. 0000e+000 


0 

3. 0024e-005 
-5. 3157e-005 
0 

-1.5239e+001 
1. 3915e-001 
-1. 4168e-002 
7 . 5402e-002 


0 

0 

0 

0 

5. 5995e~001 
0 
0 
0 


» b 


b = 


Columns 1 through 8 


-1 . 7943e-0O2 -1.7943e-002 
-8, 3378e-001 -8 . 3378e-001 
-7. 3127e+000 -7.3127e+000 
0 0 


1. 8059e-002 -1.8059e-002 
8. 3878e+000 -8.3878e+000 
-1 . 5399e-001 1.5399e-001 

0 0 


-6. 1980e-002 
-2. 9910e-001 
-2.5975e+000 
0 
0 
0 
0 
0 


-2. 2331e-002 
-5. 2179e-001 
-2. 4202e+000 
0 

4. 7930e-004 
4.3692C+000 
-2 . 3691e-D01 
0 


-2. 2331e-002 
-5. 2179e-001 
-2. 4202e+000 
0 

-4. 7930e-004 
-4. 3693e+000 
2 . 369Se-001 
0 


-1. 4240e-001 
3. 6383e-001 
2. 8081e+000 
0 

6. 1131e-002 
-4. 1643e+000 
-7. 9366e-001 
0 


-1. 4240e-001 
3. 6383e-001 
2. 8081e+000 
0 

-6, 1131e-002 
4. 1643e+000 
7. 9366e-001 
0 


-1 . 3153e-002 
1 . 8624e-002 
-1. 0550e-001 
0 

2. 3480e-002 
-1.5024e-001 
-l.Q674e-001 
0 


» c 


: ~ 




9. 9717e-001 

7. 5188e-002 

0 

0 

-4, 9201e-003 

6. 5251e-002 

0 

0 

0 

0 

l.OOOOc+OOO 

0 

0 

0 

0 

1. OOOOe+OOO 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-3. 5433e-004 

1.5489e-003 

0 

0 

3. 7342e-019 

0 

0 

0 

2. 2I30e-003 

1. 9041e-002 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6. 5436e-002 

0 

0 

0 

0 

9. 9717e-001 

7. 5188e-002 

0 

0 

-7. 5188e-002 

9. 9717e-001 

0 

0 

0 

0 

1. OOOOe+OOO 

-1.3378e-004 

0 

0 

0 

-2.S268e-004 

4. 9947e-005 

3. 7661e-006 

0 

-2.7958e-004 

-1.2376e-005 

-9.3318e-007 

0 
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d = 


Columns 1 through 8 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


-5. 5769e-004 -5.5769e-004 
5. 6129e-004 -5.6129e-004 
2. 5915e-002 2.5915e-002 


0 

0 

0 

0 

0 

0 

0 

0 

-1. 9264e-003 
0 

9. 2962e-003 


-6 . 9408e-004 -6.9408e-004 
1 . 4897e-005 -1.4897e-00S 
1 . 6218e-002 1.6218e~002 


-4. 4259e-003 -4.4259e-003 
1.9000e-003 -1.9000e-003 
■1. 1308e-002 -1.1308e-002 


0 

0 

0 

0 

0 

0 

0 

0 

-4. 0881e-004 
7. 2978e-004 
-5. 7885e-004 


Columns 9 through 11 


0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 


-4. O801e-QO4 -2 . 9148e-006 


-7. 2978e-004 0 

-5. 7885e-004 3.0557e-003 


0 

0 

0 

0 

0 

0 

0 

0 

0 

3. D557e-003 
0 
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II. A Design Challenge Solution 
A. Introduction 

One “solution” to the design challenge described in Part I of this document will 
now be presented. The controller design is similar to that presented in Ref. 1, albeit 
somewhat simpler in implementation. The design approach involves sliding mode 
control (SMC) utilizing asymptotic observers and reference model “hedging” to reduce 
the SMC system’s sensitivity to parasitic dynamics, here consisting of the actuator 
dynamics that are ignored in the design procedure. 

B. Solution 

A. Controller Description 

Figure II. 1 shows the basic control system structure. Figure II.2 is a detailed 



Figure II. 1 The SMC system structure 

representation of the continuous form of the “discretized SMC system” identified in Fig. 
II. 1 . Figure II. 3 is a more detailed representation of the “sliding mode controllers” shown 
in Fig. II. 2. Likewise, Figs. II.4, II.5 and II.6 are detailed representations of the “hedge 
system”, “reference models” and “asymptotic observers” of Fig. II. 2 It should be noted 
that all of the possible sensed quantities available with the ICE model were input to the 
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observers. The control allocation matrix was created using a pseudo-inverse design 
technique [1]. 



Figure II.2 Detail of continuous form of “discretized SMC system” in Fig. II. 1 
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Figure II.4 Detail of “hedge system” of Fig. II. 2 
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alpha edged reference model 



Figure II. 5 Detail of “reference models” of Fig. II. 2 
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observer poles 



Figure II.6 Detail of “asymptotic observers” in Fig. II.2 


Details of the design procedure will be omitted here. The reader is referred to 
Ref. 1 for a thorough discussion of the approach. It should be mentioned here, however, 
that the SMC design to be evaluated here is not an “optimum” one. For example, the 
hedging logic was created for a single flight condition (Mach No. = 0.6, Altitude = 
25,000 ft) and was left invariant for the remainder of the flight conditions evaluated. 

It should be noted that the pitch-rate command system was obtained by first 
designing an alpha-command system. Then a control stick command filter was placed in 
series with the pilot’s longitudinal stick command. The filter was a simple lag/lead 
network that produced rate-like pitch attitude responses when input to the alpha- 
command system. Since an alpha-command system was actually in evidence, no phugoid 
damping was provided by the control system. However, this produced no problems as 
the pilot control of pitch attitude provided excellent phugoid damping. 
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2.0 Solution Results 


Handling Qualities 

The bandwidth/phase delay values for each of the four flight conditions are given 

below. 

Mach No. = 0.3, Altitude = 1 5.000 ft (wings level) 



Pitch 

Roll 

“BWphase 

3.69 

2.35 

(rad/sec) 



x p 

0.071 

0.11 

(sec) 




Mach No. = 0.6. Altitude = 25,000 ft (wings level) 



Pitch 

Roll 

WBW phase 

3.69 

2.35 

(rad/sec) 



x p 

0.071 

0.11 

(sec) 




Mach No. = 0.6, Altitude = 25.000 ft (steady turn) 



Pitch 

Roll 


4.42 

2.32 

(rad/sec) 



t p 

0.065 

0.12 

(sec) 




Mach No. = 0.9. Altitude = 35,000 ft 



Pitch 

Roll 

“BWphase 

5.14 

2.32 

(rad/sec) 



Xp 

0.063 

0.12 

(sec) 
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Referring to Fig. 1.8 in Part I of this document, these values are each seen to predict Level 
1 handling qualities. 


b. Tracking Performance 

Tracking performance was examined for flight conditions for nominal and failed 
cases. In the case of failures, the failure was introduced 20 seconds into a 50 second run. 
The entire failure ensemble described in Section I.C.8 was accommodated by the 
controller. Performance was in the desired categories as described in Section I.C.7.b 
except for the following: 

Mach No. = 0.3, Altitude = 1 5,000 ft, wings-level, nominal 
only adequate performance for beta 

Mach No. = 0.6, Altitude = 25,000 ft; steady turn ,nominal 
only adequate performance for theta 

Mach No. = 0.6, Altitude = 25,000 ft; steady turn, failed 
only adequate performance for theta 


Mach No. = 0.9, Altitude = 35,000 ft, failed 

adequate performance not attained for phi 


The last of these conditions was the most serious in terms of performance, where 
adequate performance was not in evidence for the phi loop. Figures II. 7 through 11.18 
demonstrate the performance on a selected set of flight conditions and vehicle modes. 
Figure 11.18 demonstrates the significant rate limiting occurring in control effector dE3 
(left elevon) after the failure. 
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Figure II. 9 P tracking Mach No. = 0.3 Altitude = 1 5,000 ft, wings level, nominal 
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Figure II. 1 1 <j> tracking Mach No. = 0.6 Altitude = 25,000 ft, steady turn, nominal 
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60 



nominal rate limits 
before failure a 150 deg/sec 



c. Stability Margins 

The stability margin criteria of Section I.C.7.C were met in all but one case. This 
involved the yaw thrust loop in the flight condition Mach No. = 0.6, Altitude - 25,000 ft, 
steady-turn. A gain margin of only 3 dB resulted. The corresponding phase margin of 34 
degrees met the criterion. 

d. Structural Coupling Considerations 

The structural coupling criteria of Section I.C.7.d were violated in four cases. All 
three involved the magnitude of the transfer function between pilot pitch input and q*>- 
The violations were 

Mach No. = 0.3, Altitude = 15,000 ft, wings-level; -9.27 dB rather than -10 dB 

Mach No. = 0.6, Altitude = 25,000 ft, wings-level; -8.21 dB rather than -10 dB 

Mach No. = 0.9, Altitude = 35,000 ft, wings-level; -7.07 dB rather than -10 dB 

Mach No. = 0.6, Altitude = 25,000 ft, steady-turn; -8.7 dB rather that -10 dB 

Figure 11.18 shows one of the transfer functions above, with the magnitude violation 
indicated. These violations were not considered serious, as they essentially could be 
interpreted as a slight increase in the 25 rad/sec criterion frequency 



structural coupling magnitude for ro > 25 rad/sec 
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5. 


Control Activity 


Five violations of the control activity criterion of Section I.C.7.e occurred. They 
are as follows: 

Mach No. = 0.9, Altitude = 35,00 ft, level flight; commands to actuators dE2, 
dE12, dE10 and dE20 exceeded the 75% figure. 

Mach No. = 0.6, Altitude = 25,000 ft, steady turn, commands to actuator dE12 
exceeded the 75% figure. 

It should be noted that these five violations were out of a total of 11 x 4 = 44 total 
possibilities (1 1 actuator commands x 4 flight conditions). The violations were each on 
the order of 1 0% above the 75% criterion value). 

6. Off-Nominal Flight Conditions 

The off-nominal flight conditions were examined with the controller tuned for 
Mach No. = 0.6, Altitude = 25,000 ft, but with vehicle dynamics associated with Mach 
No. = 0.3, Altitude = 15,000 ft and Mach No. = 0.9, Altitude = 35,000 ft. Figures 11.19 
through 11.24 show the results. One tracking performance violation occurred for beta 
tracking for the with Mach No. = 0.3, Altitude = 15,000 ft case where only adequate 
performance could be obtained. 



Mach No. = 0.6, Altitude = 25,000 ft, but vehicle dynamics are Mach No. - 0.3, 

altitude = 15,000 ft 
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Figure 11.20 <j> tracking for off-nominal condition; controller tuned for 
Mach No. = 0.6, Altitude = 25,000 ft, but vehicle dynamics are Mach No. = 0.3 

altitude = 15,000 ft 







Figure 11.23 4> tracking for off-nominal condition; controller tuned for 
Mach No. = 0.6, Altitude = 25,000 ft, but vehicle dynamics are Mach No. = 0.9, 

altitude = 35,000 ft 
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Figure 11.24 p tracking for off-nominal condition; controller tuned for 

Mach No. = 0.6, Altitude = 25,000 ft, but vehicle dynamics are Mach No. = 0.9, 

altitude = 35,000 ft 

Scheduling the controller designed herein would be based upon Mach No. and 
Altitude. In the implementation used, the only the asymptotic observers and the control 
distribution matrix were a function of flight condition. The SMC and hedging systems 
were independent of flight condition. It should be noted that is likely that an improved 
designs could be obtained by consideration of tuning the hedging system to the flight 
condition. 

3. Discussion 

There was some concern about the amount of control activity in the design just 
presented. Therefore, for comparison, a simple controller obtained with classical loop- 
shaping techniques was designed and evaluated in the same tracking task for flight 
condition Mach No. = 0.3, Altitude = 15,000 ft. This latter design had no asymptotic 
observers, and only three variables were assumed to be measured: alphaw, ps, and beta. 
Figures 11.25 and 11.26 compare the control activity in the left elevon (dE03). As can be 
seen, the SMC design has slightly higher frequency content, but significantly less motion 
amplitude. Thus, the SMC design was considered to be acceptable from the standpoint of 
sensor noise transmission. 
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Figure 11.26 Left elevon activity (dE03) for Mach No. = 0.3, Altitude = 15,000 ft, for 

SMC design 


C. Conclusions 


The controller design that was described and exercised in Part II was intended to 
serve as a test case for the design challenge. As pointed out in this section, a number of 
violations of the design specifications occurred, although none were “show-stoppers” in 
the sense that system stability was compromised. The SMC approach allowed the 


67 



pilot/vehicle system to accommodate the entire suite of failures described in Section 
I.C.8, which is notable. 
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